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ABSTRACT 



The dramatic increase in the number of known gamma-ray pulsars since the launch 
of the Fermi Gamma-ray Space Telescope (formerly GLAST) offers the first opportu- 
nity to study a sizable population of these high-energy objects. This catalog summa- 
rizes 46 high-confidence pulsed detections using the first six months of data taken by 
the Large Area Telescope (LAT), FermVs main instrument. Sixteen previously un- 
known pulsars were discovered by searching for pulsed signals at the positions of bright 
gamma-ray sources seen with the LAT, or at the positions of objects suspected to be 
neutron stars based on observations at other wavelengths. The dimmest observed flux 
among these gamma-ray-selected pulsars is 6.0 x 10 -8 ph cm -2 s _1 (for E >100 MeV). 
Pulsed gamma-ray emission was discovered from twenty-four known pulsars by using 
ephemerides (timing solutions) derived from monitoring radio pulsars. Eight of these 
new gamma-ray pulsars are millisecond pulsars. The dimmest observed flux among the 
radio-selected pulsars is 1.4 x 10~ 8 ph cm -2 s" 1 (for E >100 MeV). The remaining six 
gamma-ray pulsars were known since the Compton Gamma Ray Observatory mission, 
or before. The limiting flux for pulse detection is non-uniform over the sky owing to 
different background levels, especially near the Galactic plane. The pulsed energy spec- 
tra can be described by a power law with an exponential cutoff, with cutoff energies in 
the range ~ 1 — 5 GeV. The rotational energy loss rate (E) of these neutron stars spans 
5 decades, from ~3 x 10 33 erg s" 1 to 5 x 10 38 erg s , and the apparent efficiencies 
for conversion to gamma-ray emission range from ~ 0.1% to ~ unity, although dis- 
tance uncertainties complicate efficiency estimates. The pulse shapes show substantial 
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diversity, but roughly 75% of the gamma-ray pulse profiles have two peaks, separated 
by > 0.2 of rotational phase. For most of the pulsars, gamma-ray emission appears to 
come mainly from the outer magnetosphere, while polar-cap emission remains plausi- 
ble for a remaining few. Spatial associations imply that many of these pulsars power 
pulsar wind nebulae. Finally, these discoveries suggest that gamma-ray-selected young 
pulsars are born at a rate comparable to that of their radio-selected cousins and that 
the birthrate of all young gamma-ray-detected pulsars is a substantial fraction of the 
expected Galactic supernova rate. 

Subject headings: catalogs - gamma rays: observations - pulsars: general - stars: neu- 
tron 

Version corrected for an Erratum sent to the Ap J, December 2010: In the original paper, an 
error was made in accounting for the delay due to interstellar dispersion in the radio phasing of 
PSR J1124— 5916. This changes the measured gamma-ray to radio lag (5) to 0.11 ± 0.01. An error 
was also made in the off-pulse phase range in Table 3 for that pulsar. This error did not affect the 
spectral results. Corrected versions of Table 3 (with the revised numbers in bold face), Figure 4, 
and Figure A20 are included here. In addition there was an error in the caption to Figure 9. The 
figure with corrected caption is included here, with the changed word in bold face. 



Introduction 



Follow i ng th e 1967 discovery of pulsars by Bell and Hewish (jHewish et al.lll968l ). iGoldl (|1968l ) 
and iPacinil (jl968l ) identified these objects as rapidly rotating neutron stars whose observable emis- 
sion is powered by the slow-down of the rotation. With their strong electric, magnetic, and gravita- 
tional fields, pulsars offer an opportunity to study physics under extreme conditions. As endpoints 
of stellar evolution, these neutron stars, together with their associated supernova remnants (SNRs) 
and pulsar wind nebulae (PWNe), help probe the life cycles of stars. 



Over 1800 rotation-powered pulsars are now listed in the ATNF pulsar catalog ({Manchester et al 

as illustrated in Figure [TJ The vast majority of these pulsars were discovered by radio tele- 
scopes. Small numbers of p ulsars have also been seen in the optical band, with more in the X-ray 
bands (see e.g. |Beckerj |2009j) . 



In the high-energy gamma-ray domain (> 30 MeV) the first ind ications for pulsar em ission 

), and 



Browning et al 



1971 



were obtained for the Crab pul sar by balloon-born e detectors (e.g. 

confirmed b y the SAS-2 satellite (|Kniffen et al.lll974h . which also found gamma radiation from the 
Vela pulsar ({Thompson et al.lll975l ). The COS-B satellite provided additional details about these 



1 http://www.atnf. csiro.au/research/pulsar/psrcat 
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two gamma-ray pulsars, including a confirmation that the Vela pulsar gamma-ray emission was not 
in phase with the radio nor d id it have the same emission pattern (light curve) as seen in the radio 
(see e.g. iKanbach et alJll98Q ). 



The Compton Gamma Ray Observatory (CGRO) expanded the number of gamma-ray pulsars 
to at least seven, with six clearly seen by the CGRO high-energy instrument, EGRET. This gamma- 
ray pulsar population allowed a search for trends, such as the increase i n efficiency rj = L^/E with 
decreasing values of the open field line voltage of the pulsar, first noted bv lAronsI (|1996l ). for gamma- 
ray luminosity L 7 and spin-down lumin osity E. A summary of gamma-ray pulsar results in the 
CGRO era is given by iThompsonl (|20Q4j ) . 



The third EGRET catalog (3EG; lHartman et al.lll999l ) included 271 sources of which ~170 re- 
mained unidentified. Determining the nature of these unidentified sources is one of the outstanding 
problems in high-energy astrophysics. Many of them are at high Galactic latitude and are most 
likely active galactic nuclei or blazars. However, most of the sources at low Galactic latitudes (\b\ < 
5°) are associated with star- forming regions and hence may be pulsars, PWNe, SNRs, winds from 
mass i ve stars, or high-ma ss X-ray binaries (e.g. iKaaret fc Cottaml Il996l ; lYadigaroglu fc Romani 
19971 ; iRomero et al.lll999l). A numb er of radio pulsars were subsequently discovered in EGRET 
error boxes (e.g. iKramer et al.ll2003l ). but gamma-ray pulsations in the archival EGRET data were 
never clearly seen. Solving the puzzle of the unidentified so urces will constra i n puls a r emission mod- 



els: p ulsar population synthesis studies, s uch as those by ICheng &: Zhangl (| 19981 ) , iGonthier et al 



( 2002 ). and McLaughlin Sz Cordes (j200ol ). indicate that the number of detectable pulsars in ei- 
ther EGRET or Fer mi data, as well as the expected ratio of radio-loud and radio-quiet pulsars 
(jHarding et al.ll2007l ). strongly depends on the assumed emission model. 



The Large Area Telescope (LAT) on the Fermi Gamma-ray Space Telescope has provided a 
major increa se in the known ga mma-ray pulsar population, inclu ding pulsars discove red first in 
gamma-rays (jAbdo et al.l l2009d ) and millisecond pulsars (MSPs) (jAbdo et al.1 l2009bl ). The first 
aim of this paper is to summarize the properties of the gamma-ray pulsars detected by Fermi-LAT 
during its first six months of data taking. The second primary goal is to use this gamma-ray pulsar 
catalog to address astrophysical questions such as: 



1 . Are all the gamma-ray pulsars consistent with one type of emission model? 

2. How do the gamma-ray pulsars compare to the radio pulsars in terms of physical properties 
such as age, magnetic field, spin-down luminosity, and other parameters? 

3. Are the trends suggested by the CGRO pulsars confirmed by the LAT gamma-ray pulsars? 

4. Which of the LAT pulsars are associated with SNRs, PWNe, unidentified EGRET sources, 
or TeV sources? 



The structure of this paper is as follows: Section 2 describes the LAT and the pulsar data 
analysis procedures; Section 3 presents the catalog and derives some population statistics from our 
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sample; Section 4 studies the LAT sensitivity for gamma-ray pulsar detection, while in Section 5 
the implications of our results are briefly discussed. Finally, our conclusions are summarized in 
Section 6. 



2. Observations and Analysis 



The Fermi Gamma-ray Space Telescope was successfully launched on 2008 June 11, carry- 
ing two gamma-ray instruments: the LAT and the Gamma-ray Burst M onitor (GBM). The LAT, 
FermVs main instru ment, is desc ri bed in detail in lAtwood et al.l (|2009l ). with early on-orbit per- 
formance reported in lAbdo et al.l (|2009ol ). It is a pair-production telescope composed of a 4 x 4 
grid of towers. Each tower consists of a silicon-strip detector and a tungsten- foil tracker /converter, 
mated with a hodoscopic cesium-iodide calorimeter. This grid of towers is covered by a segmented 
plastic scintillator anti-coincidence detector. The LAT is sensitive to gamma rays with energies in 
the range from 20 MeV to greater than 300 GeV, and its on-axis effective area is ~ 8000 cm 2 for 
E > 1 GeV. The gamma-ray point spread function (PSF) is energy dependent, and 68% of photons 
have reconstructed directions within 6q% ~ 0°8E~°' S of a point source, with E in GeV, leveling off 
to #68 ^ 0?1 for E > 10 GeV . Effective area, PSF, and energy resolution are tabulated into bins of 
photon energy and angle of incidence relative to the LAT axi s . The tables are called "instrument 
response functions, and are described in detail in lAbdo et al.1 (|2009ol ). This work uses the version 
called P6_v3_dif f use. 

Gamma-ray events recorded with the LAT have time stamps that are derived from a GPS- 
synchroniz ed clock on board t he Fermi satellite. The accuracy of the time stamps relative to UTC 
is < 1 fis (jAbdo et al.ll2009ol ). The timing chain from the GPS-based satellite clock through the 
barycentering and epoch folding software has been shown to be accurate to be tter than a few /is 
for binary orbits, and significantly better for isolated pulsars (jSmith et al.l 120081 ). 



The LAT field-of-view is about 2.4 sr. Nearly the entire first year in orbit has been dedi- 
cated to an all-sky survey, imaging the entire sky every two orbits, i.e. every 3 hours. Any given 
point on the sky is observed roughly l/6 th of the time. The LAT's large effective area and ex- 
cellent source localization coupled with improved cosmic-ray rejection led to the detection of 46 



gamma-ray pulsars in the first six mo nths of LA^ 
ray pulsa rs clearly seen with EGRET ( Thompson! 



EGRET (IRamanamurthy et al.lll996l : iKaspi et al 



obs ervations. These include the six gamma- 



2004 ). two young pulsars seen marginally with 



2000), the MSP seen mar ginally with EGRET 



(|Kuiper et alihoooh . PSR J2021+3651 discovered in ga mma-rays by AGILE (jHalpern et al.ll2008h 
and some of the other pulsars also studied by AGILE (jPellizzoni et al.ll2009l ) 



During the LAT commissioning period, several configuration settings were tested that affected 
the LAT energy resolution and reconstruction. However, these changes had no effect on the LAT 
timing. Therefore, for the spectral analyses, the data were collected from the start of the Fermi 
sky-survey observations (2008 August 4, shortly before the end of the commissioning period) until 
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2009 February 1, while the light curve and periodicity test analysis starts from the first events 
recorded by the LAT after launch (2008 June 25) and also extends through 2009 February 1. 



2.1. Timing Analysis 

We have conducted two distinct pulsation searches of Fermi LAT data. One search uses the 
ephemerides of known pulsars, obtained from radio and X-ray observations. The other method 
searches for periodicity in the arrival times of gamma-rays coming from the direction of neutron 
star candidates ("blind period searches"). Both search strategies have advantages. The former is 
sensitive to lower gamma-ray fluxes, and the comparison of phase- aligned pulse profiles at different 
wavelengths is a powerful diagnostic of beam geometry. The blind period search allows for the 
discovery of new pulsars with selection biases different from those of radio searches, such as, for 
example, favoring pulsars with a broader range of inclinations between the rotation and magnetic 
axes. 

For each gamma-ray event (index i), the topocentric gamma-ray arrival time recorded by the 
LAT is transferred to times at the solar-system barycenter ti by correcting for the position of Fermi 
in the solar-system frame. The rotation phase (fii(ti) of the neutron star is calculated from a timing 
model, such as a truncated Taylor series expansion, 

Uu) = <Po + }_^ — 7-— rv — ■ w 

j=0 VJ -r J- 

Here, To is the reference epoch of the pulsar ephemeris and 4>q is the pulsar phase at t = Tq. The 
coefficients fj are the rotation frequency deri yatives of order j. Th e rotation period is P = l//o- 



Different timing models, described in detail in lEdwards et al.l (|2006l ). can take into account various 



physical effects. Most germane to the present work is accurate <f>i(ti) computations, even in the 
presence of the rotational instabilities of the neutron star called "timing noise" . "Phase- folding" a 
light curve, or pulse profile, means filling a histogram with the fractional part of the 0, values. An 
ephemeris includes the pulsar coordinates necessary for barycentering, the fj and Tq values, and 
may include parameters describing the pulsar proper motion, glitch epochs, and more. The radio 
dispersion measure (DM) is used to extrapolate the radio pulse arrival time to infinite frequency, 
and uncertainties in the DM translate to an uncertainty in the phase offset between the radio and 
gamma-ray peaks. 



2.1.1. Pulsars with Known Rotation Ephemerides 



The ATNF database (version 1.36. iManchester et al.ll2005l ) lists 1826 pulsars, and more have 
been discovered and await publication (Figure [I]) . The LAT observes them continuously during 
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the all-sky survey. Phase-folding the gamma rays coming from the positions of all of these pulsars 
(consistent with the energy-dependent LAT PSF) requires only modest computational resources. 
However, the best candidates for gamma-ray emission are the pulsars with high E, which often 
have substantial timing noise. Ephemerides accurate enough to allow phase-folding into, at a 
minimum, 25-bin phase histograms can degrade within days to months. The challenge is to have 
contemporaneous ephemerides. 

We have obtained 762 contemporaneous pulsar ephemerides from radio observatories, and 5 
from X-ray telescopes, in two distinct groups. The first consists of 218 pulsars with high spin-down 
power (E > 10 34 erg s" 1 ) timed regu larly as part of a c ampaign by a consortium of astronomers 



for the Fermi mission, as described in lSmith et al.l (|2008l ). With one exception (PSR J1124— 5916, 
which is a faint radio source with large timing noise), all of the 218 targets of the campaign have 
been mo nitored since shortly bef ore Fermi launch. Some results from the timing campaign can be 
found in IWeltevrede et al.l (|2009l ). 



The second group is a sample of 544 pulsars from nearly the entire P — P plane (Figure [2]) 
being timed for other purposes for which ephemerides were shared with the LAT team. These 
pulsars reduce possible bias of the LAT pulsar searches created by our focus on the high E sample 
requiring frequent monitoring. Gamma-ray pulsations from six radio pulsars with E < 10 erg s _1 
were discovered in this manner, all of which are MSPs. 

Table [1] lists properties of the 46 detected gamma-ray pulsars. Five of the 46 pulsars, all MSPs, 
are in binary sys tems. The period first derivative P is corrected for the kinematic Shklovskii effect 
(|Shklovskiilll970h : P = P ohs - /i 2 P ohs d/c, where [i is the pulsar proper motion, and d the distance. 
The correction is small except for a few MSPs (jAbdo et al.l l2009bl ) . The characteristic age is 
r c = P/2P and the spin-down luminosity is 



E = 4nt*IPP~ 



(2) 



taking the neutron star's moment of inertia I to be 10 45 g cm 2 . The magnetic field at the light 
cylinder (radius i?Lc = cP/2ir) is 



3I8tt 4 P 
c 3 P 5 



1/2 



2.94 x 10 8 (PP" 5 ) 1/2 G. 



(3) 



Table [2] lists which observatories provided ephemerides for the gamma-ray pulsars. An "L" 
indicates that the pulsar was timed using LAT gamma rays, as described in the next Section. 



"P" is the Parkes radio telescope ([Manchester! 120081 ; IWeltevrede et al.ll2009l ). The majority of 
the Parkes observations were carried out at intervals of 4 - 6 weeks using the 20 cm Multibeam 



receiver (jStavelev-Smith et al.lll996l ) with a 256 MHz band centered at 1369 MHz. At 6 month 
intervals observations were also made at frequencies near 0.7 and 3.1 GHz with bandwidths of 32 
and 1024 MHz respectively. The required frequency resolution to avoid dispersive smearing across 
the band was provided by a digital spectrometer system. 
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"N" is the Nangay radio telescope (jTheureau et al.ll2005l ). Nancay observations are carried 
out every three weeks on average. The recent version of the BON backend is a GPU-based coher- 
ent dedispersor allow ing the processing of a 128 MHz bandwidth over two complex polarizations 
(jCognard et al.l 12009). The majority of the data are collected at 1398 MHz, while for MSPs in 
particular, observations are duplicated at 2048 MHz to allow DM monitoring. 

"J" is th e 76-m Lovell radio telescope at Jodrell Bank in the United Kingdom. Jodrell Bank 
observations (jHobbs et al.l [2004J) were carried out at typical intervals of between 2 days and 10 
days in a 64- MHz band centered on 1404 MHz, using an analog filterbank to provide the frequency 
resolution required to remove interstellar dispersive broadening. Occasionally, observations are also 
carried out in a band centered on 610 MHz to monitor the interstellar dispersion delay. 

"G" is the 100-m NRAO Green Bank Telescope (GBT). PSR J18 33-1034 was observ ed monthly 
at 0.8 GHz with a bandwidth of 48 MHz using the BCPM filter bank ^Backer et al.lll997h . The other 
pulsars monitored at the GBT were observed every two weeks at a frequency of 2 GHz across a 
600 MHz band with the Spigot spectrometer (jKaplan et al.l 120051 ) . Individual integration times 
ranged between 5 minutes and 1 hour. 

Arecibo ("A") observations of the very faint PSRs J1930+1852 and J2021+3651 are carried 
out every two weeks, with the L-wide receiver (1100 to 1730 MHz). The back-ends used are the 
Wideband Arecibo Pulsar Processor (WAPP) correlators (Dowd et al. 2000), each with a 100 
MHz-wide band. The antenna voltages are 3-level digitized and then auto-correlated with a total 
of 512 lags, accumulated every 128 fis, and written to disk as 16-bit sums. Processing includes 
Fourier-transforms to obtain power spectra, which are then dedispersed and phase-folded. The 
average pulse profiles are cross-correlated with a low-noise template profile to obtain topocentric 
times of arrival. 

"W" is the Westerbork Synthesis Radio Telescope, with which observations were made approx- 
imately monthly at central frequencies of 328, 382 and 1380 MHz with bandwidths of 1 MHz at the 



lower frequencies and 80 MHz at the higher frequencies. The PuMa pulsar backend (jVoute et al 



20021 ) was used to record all the observations. Folding and dedispersion were performed offline. 



The rms of the radio timing residuals for most of the solutions used in this paper is < 0.5% 
of a rotation period, but ranges as high as 1.2% for five pulsars. This is adequate for the 50- or 
25-bin phase histograms used in this paper. The ephemerides used for this catalog will be available 
on the Fermi Science Support Center data serverio 



2 http://fermi. gsfc.nasa.gov/ssc/ 
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2.1.2. Pulsars Discovered in Blind Periodicity Searches 



For all 16 of the pulsars found in the blind searches of the LAT data, we determined the timing 
ephemerides used in this catalog directly from the LAT data as described below. In addition, for 
two other pulsars the LAT data provided the best available timing model. The first is the radio- 
quiet pulsar Geminga. Since Geminga is such a bright gamma-ray pulsar, it is best timed directly 
using gamma-ray observations. During the period b etween EGRET and Ferm i, occasional XMM- 
Newton observations maintained the timing model (j Jackson &: Halpernl 120051) but a substantially 
improved ephemeris has now been derived from the LAT data (lAbdo et al.ll2009al ). The second is 
PSR J1124— 591 6, which is extremel y faint in the radio (see Table [1]) and exhibits a large amount 
of timing noise (jCamilo et al.ll2002d ). In this Section, we briefly describe the blind pulsar searches 
and how the timing models for these pulsars are created. These pulsars have an "L" in the "ObsID" 
column of Table [2j 

Even though L 7 of a young pulsar can be several percent of E, the gamma-ray counting rates 
are low. As an example, the LAT detects a gamma ray from the Crab pulsar approximately every 
500 rotations, when the Crab is well within the LAT's field-of-view. Such sparse photon arrivals 
make periodicity sear c hes difficult. Exten sive searches for pulsations performed on EGRET data 
([Chandler et al.l l200ll ; IZiegler et al.l 120081 ) were just sensitive enough to detect the very bright 
pulsars Vela, Crab, and Geminga in a blind search, had they not already been known as pulsars. 
Blind periodicity searches of all other EGRET sources proved fruitless. 

By contrast, the improvements afforded by the LAT have enabled highly successful blind 
searches for pulsars. In the first six months of operati on, we discovered a total of 16 new pulsars 
in direct pulsation searches of the LAT data (see e.g. Abdo et al. 2008 . 2009c!) . A computation - 
ally efficient time-difference search technique made these searches possible (jAtwood et al.ll2006l ). 
enabling searches of hundreds of Fermi sources to be performed on a small computer cluster with 
only a modest loss in sensitivity compared to fully coherent search techniques. Still, owing to the 
large number of frequency and frequency derivative trials required to search a broad parameter 
space, the minimum gamma-ray flux needed for a statistically significant detection is considerably 
higher than the minimum flux needed for the phase-folding technique using a known ephemeris (as 
in Section 2.1, Eq. 1). 

We performed these blind searches on ~ 100 candidate sources identified before launch and on 
another ~ 200 newly detected LAT sources. The parameter space covered by the blind searches 
included frequencies from 0.5 Hz to 64 Hz (periods of 156.25 ms to 2 s), and a frequency derivative 
from zero to the spin-down of the young Crab pulsar (f-\ = -3.7 x 10 -10 ), which covers ~ 86% of 
the pulsars contained in the ATNF database (jAbdo et al.ll2009d ). Of the 16 pulsars detected in 
these searches, 13 are associated with previously known EGRET sources. The discoveries include 
several long-suspected pulsars in SNRs and PWNe. 



These 16 pulsars are gamma-ray selected, as they were discovered by the LAT and thus the 
population is subject to very different selection effects than the general radio pulsar population. 
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However, this does not necessarily imply that they are radio quiet. For several cases, deep radio 
searches have already been performed on known PWNe or X-ray point sources suspected of har- 
boring pulsars. In most cases, new radios searches are required to investigate whether there is a 
radio pulsar counterpart down to a meanin gful luminosity limi t. These searches are now being 
undertaken and are yielding the first results (|Camilo et al.ll2009l ). 



For these 18 pulsars (16 new plus Geminga and PSR J1124— 5916), we derived timing models 
from the LAT data using the procedur e summarized he re. A more detailed description of pulsar 
timing using LAT data can be found in iRav et al.l (|2010l ) . 



The LAT timing analysis starts from the first events recorded by the LAT after launch (2008 
June 25) and extends through about 2009 May 1. During the commissioning period, several config- 
uration settings were tested that affected the LAT energy resolution and reconstruction. However, 
these changes had no effect on the LAT timing. We selected photons from a small radius of interest 
(ROI) around the pulsar of < 0.5° or < 1° (see further Section 2.1.3 and Table [2|). For this pulsar 
timing analysis, we used diffuse class photons with energies above a cutoff (typically E > 300 MeV) 
selected to optimize the signal-to-noise ratio for that particular pulsar. We converted the photon 
arrival times to the geocenter using the gtbary science tool 2 . This correction removes the effects 
of the spacecraft motion about the Earth, resulting in times as would be observed by a virtual 
observatory at the geocenter. 



Using an initial timing model for the pulsar, we then used Tempo2 (jHobbs et al.ll2006l ) in its 
predictive mode to generate polynomial coefficients describing the pulse phase as a function of time 
for an observatory at the geocenter. Using these predicted phases, we produced folded pulse profiles 
over segments of the LAT observation. The length of the segments depends on the brightness of the 
pulsar but are typically 10-20 days. We then produced a pulse time o f arrival (TO A) for each data 
segment by Fourier domain cross-correlation with a template profile (lTaylorlll992l ). The template 
profile for most of the pulsars is based on a multi-Gaussian fit to the observed LAT pulse profile. 
However, in the case of Geminga, which has very high signal-to- noise and a complex profile not 
well described by a small number of Gaussians, we used a template profile that was the full mission 
light curve itself. 

Finally, we used Tempo2 to fit a timing model to each pulsar. For most of the pulsars, the 
model includes pulsar celestial coordinates, frequency and frequency derivative. In several cases, 
the fit also required a frequency second derivative term to account for tim ing noise. In the case 
of PSR J1124— 5916, we required three sinusoidal terms (jHobbs et al.ll2006l ) to model the effects of 
the strong timing noise in this source. Two pulsars (J1741— 2054 and J1809— 2332) have positions 
too close to the ecliptic plane for the Declination to be well-constrained by pulsar timing and thus 
we fixed the positions based on X-ray observations of the presumed counterpart^. For Geminga 



3 Pulsar timing positions are measured by fitting the sinusoidal delays of the pulse arrival times associated with the 
Earth moving along its orbit. For pulsars very close to the ecliptic plane the derivative of this delay with respect to 
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and PSR J1124— 591 6 we also used known, fixed positions (jCaraveo et al,lll998l ; lFahertv et al.ll2007l : 



Camilo et al.ll2002d j because they were of much higher precision than could be determined from 
less than one year of Fermi timing. The rms of the timing residuals are between 0.5 and 2.9% of a 
rotation period, with the highest being for PSR J1459— 60. The Tempo2 timing models used for 
the catalog analysis will be made available online at the FSSC web site 2 . 



2.1.3. Light Curves 

The light curves of 46 gamma-ray pulsars detected by the LAT are appended to the end of 
this paper, in Figures lA-ll to lA-461 The gray light curve in the top panel includes all photons with 
E > 0.1 GeV, while the other panels show the profiles in exclusive energy ranges: E > 1.0 GeV 
(with E > 3.0 GeV in black) in the second panel from the top; 0.3 to 1.0 GeV in the next panel; 
and 0.1 to 0.3 GeV in the fourth panel. Phase-aligned radio profiles for the radio-selected pulsars 
are in the bottom panel. The light curves are plotted with N = 25 or 50 bins, with 25 bins used 
when required to keep at least 50 counts per bin in the peak of the light curve or to prevent undue 
smearing due to the accuracy of the timing model. 

Table [3] lists light curve parameters, taken from the > 100 MeV profiles (top panel of Figures 
A-l-A-46). For some pulsars (e.g. PSR J1420— 6048) the peak multiplicity is unclear: the data are 
consistent with both a single, broad peak and with two closely spaced narrow peaks. For the Crab 
(PSR J0534+2200) the phase offset of the first gamma-ray peak is relative to radio "precursor" 
peak, not visible in Figure A-8. 



Table [2] lists the Z\ (IBuccheri et al.lll983l ) and H (jde Jager et al.lll989l ) periodicity test values 
for E > 0.3 GeV. Only one trial is made for each pulsar, and the significance calculations do not take 
into account the trials factor for the ~ 800 pulsars searched. Detection of gamma-ray pulsations 
are claimed when the significance of the periodicity test exceeds 5a (i.e. a chance probability of 
< 6 x 10~ 7 ). We have used the Z-test with m = 2 harmonics (Z$) which provides an analytical 
distribution function for the null hypothesis described by a x 2 distribution with 2m degrees of 
freedom. The H-test uses Monte-Carlo simulations to calculate probabilities, limited to a minimum 
of 4 x 10~ 8 (equivalent to 5.37a). Each method is sensitive to different pulse profile shapes. Four 
pulsars in the catalog fall short of the 5<7 significance threshold in the six-month data set with the 
selection cu t s appl ied here: the 3 MSPs J0218+4232, J0751+1807, and J1744-1134 reported in 



Abdo et al.l (j2009bl ) . and the radio pulsar PSR J2043+2740. The characteristic pulse shape as well 



as the trend of the significance versus time lead us to include these four in the catalog. 

Table [2] also lists "maxROI" , the maximum angular radius around the pulsar position within 
which gamma-ray events were kept, generally 1?0, but 0?5 in some cases. The choice was made 



ec liptic latitude is g reatly reduced and thus such pulsars have one spatial dimension poorly constrained, as discussed 
in lRav et all l|201oh . 
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by using the energy spectrum for the phase-averaged source, described in Section 2.2, to maximize 
S 2 /N over a grid of maximum radii and minimum energy thresholds (where S is the number of 
counts attributed to the point source, and N is the number of counts due to the diffuse background 
and neighboring sources). We selected photons within a radius #68 of the pulsar position, requiring 
a radius of at least 0?35, but no larger than the reported maxROI. 

We estimated the background level represented by the dashed horizontal lines in Figures IA-11 
to lA-46l from an annulus between 1°< 6 <2° surrounding the source. Nearby sources were removed, 
and we normalized to the same solid angle as the source ROI. The poor spatial resolution of the 
LAT at low energies can blur structured diffuse emission and bias this background estimate. The 
levels shown are intended only to guide the eye. Detailed analyses of off-pulse emission will be 
discussed in future work. 



2.2. Spectral Analysis 



The pulsar spectra were fitted with an exponentially cutoff power-law model of the form 



dN 
dE 



KE GeV exp 



E 



^cutoff 



(4) 



in which the three parameters are the photon index at low energy T, the cutoff energy -^cutoff) an d 
a normalization factor K (in units of ph cm " 2 s _1 MeV -1 ), in keeping with the observed spectral 



shape of bright pulsars (jAbdo et al.l l2009el ) . The energy at which the normalization factor K is 
defined is arbitrary. We chose 1 GeV because it is, for most pulsars, close to the energy at which 
the relative uncertainty on the differential flux is minimal. 

We wish to extract the spectra down to 100 MeV in order to constrain the power-law part of 
the spectrum, and to measure the flux above 100 MeV directly. Because the spatial resolution of 
Fermi is not very good at low energies (~ 5° at 100 MeV), we need to account for all neighboring 
sources and the diffuse emission toge ther with each pu lsar. This was done using the framework 
used for the LAT Bright Source List (jAbdo et al.l 12009 Jl. A 6-month s ource list was generated in 
the same way as the 3-month source list described in lAbdo et al. ( 2009i ). but co vering the extende d 
period of time used for the pulsar analysis. We have added the source Cyg X-3 (jAbdo et al.ll2009ll ). 
although it was not detected automatically as a separate source, because it is very close to PSR 
J2032+4127, and impacts the spectral fit of the pulsar. Cyg X-3 was fit with a simple power-law 
as were all other non-pulsar sources in the list. 



We used a Galactic diffuse model desig nated 54_77Xvarh7S calculated using GALPROPjj, an 



evolution of that used in lAbdo et all (j2009il ). A similar model, gll_iem_v02, is publicly available 2 . 



We kept events with E > 100 MeV belonging to the diffuse event class, which has the tightest 



4 http: / /galprop. stanford.edu/ 
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cosmic-ray background rejection (jAtwood et al.l 12009 ) . To avoid contamination by gamma-rays 



produced by cosmic-ray interactions in the Earth's atmosphere, we select time intervals when the 
entire ROI, of radius 10° around the source, has a zenith angle < 105°. We extracted events in a 
circle of radius 10° around each pulsar, and included all sources up to 17° into the model (sources 
outside the extraction region can contribute at low energy). Sources further away than 3° from 
the pulsar were assigned fixed spectra, taken from the all-sky analysis. Spectral parameters for the 
pulsar and sources within 3° of it were left free for the analysis. 

The fit was performed b y maximizing unbin ned likelihood (direction and energy of each event 
is considered) as described in Abdo et al. ( 2009i ) and using the minuit fitting engin^j. The uncer- 



tainties on the parameters were estimated from the quadratic development of the log(likelihood) 
surface around the best fit. In addition to the index T and the cutoff energy E cuto ff which are 
explicit parameters of the fit, the important physical quantities are the photon flux iqoo (in units 
of ph cm -2 s _1 ) and the energy flux Gioo (in units of erg cm -2 s _1 ), 

,100 Gcv dN 

F wo = / — dE, and (5) 

J 100 MeV atj 

,100 GeV dN 

G wo = / E—dE. (6) 

JlOOMeV atj 

These derived quantities are obtained from the primary fit parameters. Their statistical uncertain- 
ties are obtained using their derivatives with respect to the primary parameters and the covariance 
matrix obtained from the fitting process. 

For a number of pulsars, an exponentially cutoff power-law spectral model is not significantly 
better than a simple power-law. We identified these by computing TS CVL toS = 2Alog(likelihood) 
(comparable to a x 2 distribution with one degree of freedom) between the models with and without 
the cutoff. Pulsars with TS'cutofr < 10 have poorly measured cutoff energies. T5 cuto ff is reported in 
Table H 

The above analysis yields a fit to the overall spectrum, including both the pulsar and any 
unpulsed emission, such as from a PWN. To do better we split the data into on-pulse and off-pulse 
samples and modeled the off-pulse spectrum by a simple power-law. The off-pulse window used for 
this background estimation is defined in the last column of Table [3j 

In a second step we re-fitted the on-pulse emission to the exponentially cutoff power-law as 
before, with the off-pulse emission (scaled to the on-pulse phase interval) added to the model and 
fixed to the off-pulse result. In many cases the off-pulse emission was not significant at the 5er or 
even 3a level, but we kept the formal best fit anyway, in order to not bias the pulsed emission 
upwards. The results summarized in Tabled] come from this on-pulse analysis. 

Using an off-pulse pure power law is not ideal for the Crab or any other PWN with synchrotron 



5 http://lcgapp.cern.ch/project/cls/work-packages/mathlibs/minuit /doc/doc. html 
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and inverse Compton components within the Fermi energy range. Judging from the Crab pulsar, 
using a simple power-law to model the off-pulse emission mainly affects the value of the cutoff 
energy. The analysis specific to the Crab, with a model adapted to the pulsar synchro tron compo 



nent lo w energies and to the high energy nebular component, yields -^cutoff ~ 6 GeV (jAbdo et al 



2010dl ). This is the value listed in Table HI The cutoff value obtained with the simplified model 
applied to most pulsars in this paper is higher (> 10 GeV). The photon and energy fluxes given 
by the two analyses are within 10% of each other. Additional exceptions in Table |4] are for PSRs 
J1836+5925 and J2021+4026. The off-pulse phase definition for these pulsars is unclear, so the 
spectral parameters reported in the Table are from the initial, phase- averaged spectral analysis. 

We have checked whether our imperfect knowledge of the Galactic diffuse emission may impact 
the pulsar paramet ers by applying the same analysis with a different diffuse model, as was done in 



Abdo et al .1 (|2009il ). The phase- averaged emission is affected. Seven (relatively faint) pulsars see 
their flux move up or down by more than a factor of 1.5. On the other hand, the pulsed flux is 
much more robust, because the off-pulse component absorbs part of the background difference, and 
the source-to-background ratio is better after on-pulse phase selection. Only two pulsars see their 
pulsed flux move up or down by more than a factor of 1.2, and none shift by more than a factor of 1.4 
when changing the diffuse model. Overall, the systematic uncertainties due to the diffuse model on 
the fluxes Fioo and Gioo, on the photon index T, and on -^cutoff scale with the statistical uncertainty. 
Adding the statistical and systematic errors in quadrature amounts, to a good approximation, to 
multiplying the statistical errors on iqoo> Gioq, and T by 1.2. The uncertainties listed in Tabled] 
and plotted in the figures include this correction. The increase in the uncertainty on E cnto s due to 
the diffuse model is < 5% and is neglected. 

Systematic uncertainties on the LAT effective area are of order 5% near 1 GeV, 10% below 0.1 
GeV, and 20% above 10 GeV. To propagate their effect on the spectral parameters in Equation 4, we 
modify the instrument response functions to bracket the nominal values, and repeat the likelihood 
calculations. This is report ed in detail for mos t of the individual LAT pulsars already referenced. 



The bias values reported in lAbdo et al.1 (|2009bl ) well describe our current knowledge of the effect of 
the uncertainties in the instrument response functions on the spectral parameters: they are ST = 
(+0.3, -0.1), 5E cutoS = (+20%, -10%), 5F W0 = (+30%, -10%), and 5G W0 = (+20%, -10%). 
The bias on the integral energy flux is somewhat less than that of the integral photon flux, due to 
the weighting by photons in the energy range where the effective area uncertainties are smallest. 
We do not sum these uncertainties in quadrature with the others, since a change in instrument 
response will tend to shift all spectral parameters similarly. 



The pulsar spectra were also evaluated using an unfolding method (|D'Agostinilll995l : iMazziotta 



20091 ). that takes into account the energy dispersion introduced by the instrument response function 
and does not assume any model for the spectral shapes. "Unfolding" is essentially a deconvolution 
of the observed data from the instrument response functions. For each pulsar we selected photons 
within 68% of the PSF with a minimum radius of 0?35 and a maximum of 5°. 
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The observed pulsed spectrum was built by selecting the events in the on-pulse phase interval 
and subtracting the events in the off-pulse interval, properly scaled for the phase ratio. The 
instrument response function, expressed as a smearing matrix, was evaluate d using the LAT Geant^- 
basecjf] Monte Carlo simulation package called Gleam ( Boinee et al, 2003 ). taking into account the 
pointing history of the source. 

The true pulsar energy spectra were th en reconstructed from the observed ones using an 
iterative procedure based on Bayes' theorem (jMazziottal 12009 ) . Typically, convergence is reached 
after a few iterations. When the procedure has converged, both statistical and systematic errors 
on the observed energy distribution can be easily propagated to the unfolded spectra. The results 
obtained from the unfolding analysis were consistent within errors with the likelihood analysis 
results. 



3. The LAT Pulsar Sample 

We describe here the astronomical context of the observed LAT pulsars, including our current 
best understanding of the source distances, the Galactic distribution and possible associations. 
We also note correlations among some observables which may help probe the origin of the pulsar 
emission. 



3.1. Distances 

Converting measured pulsar fluxes to radiated power requires reliable distance estimates. An- 
nual trigonometric parallax measurements are the most reliable, but are generally only available 
for a few relatively nearby pulsars. 

The most commonly used technique to obtain radio pulsar distances exploits the pulse de- 
lay as a function of wavelength by free electrons along the path to Earth. A distance can be 
compu ted from the DM cou pled to an electron density distribution model. We use the NE2001 



model (jCordes fc Lazioll2002l ) unless noted otherwise. It assumes uniform electron densities in and 
between the Galactic spiral arms, with smooth transitions between zones, and spheres of greater 
density for specific regions such as the Gum nebula, or surrounding Vela. Specific lines-of-sight can 
traverse unmodeled regions of over- (or under-) density, as, for example, along the tangents of the 
spiral arms, causing significant discrepancies between the true pulsar distances and those inferred 
from the electron-column density. 

A third method, kinematic, associates the pulsar with objects whose distance can be measured 
from the Doppler shift of absorption or emission lines in the neutral hydrogen (HI) spectrum, 



http://geant4.web.cern.ch/geant4/ 
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together with a rotation curve of the Galaxy. It breaks down where the velocity gradients are 
very small or where the distance-velocity relation has double values. The associations are often 
uncertain, and these distance measurements can be controversial. 

In a small number of cases, the distance is evaluated either from X-ray measurements of the 
absorbing column at low energies (below 1 keV), or from consideration of the X-ray flux assuming 
some standard parameters for the neutron star. 

Table [5] presents the best known distances of 37 pulsars detected by Fermi, the methods used 
to obtain them, and the references. For distances obtained from the NE2001 model and the DM, 
the reference indicates the DM measurement. We assume a minimum DM distance uncertainty 
of 30%. When distances from different methods disagree and no method is more convincing than 
the other, a range is given, and 30% uncertainties on the upper and lower values are used. For 
the remaining 9 Fermz-discovered pulsars no distance estimates have been established so far. Here 
follow comments for some of the distance values reported in Table [5j 



PSR J0205+6449 - The pulsar is in the PWN 3C 58. NE2001 gives 4.5 kpc for DM=141 



cm 



band reddening (jFesen et al 



1988 



~ 3 pc in this direction ( Camilo et al. 2002d 'l. Using HI ab sorption and emissio n lines from the 
PWN yields fro m 2.6 kpc dGreen fc Gulll Il982h to 3.2 kpc rtRoberts et al.l ll993h. The l ower V- 



20081 ) compared to the Galactic-disk edge (jSchlegel et al.lll998l ) 
suggests that the P WN is in the rang e 3-4 kpc. Table [5] quotes the distance range found by 



Green & Gulll (| 19821 ) and iRoberts et al.1 (|1993T ) 



PSR J0218+4232 - The DM measurements from iNavarro et all (j 19951 ) together with NE2001 
yield 2.7 ± 0.8 kpc. Comparing the pulsar characteristic age with the coo ling models of its white- 
dwarf companion gives a distance range of 2.5 to 4 kpc (IBassa et al.ll2003l ). 



PSR J0248+ 6021 - The DM of 376 cm" 3 pc (ICoenard et al.ll2010h puts this pulsar beyond 
the edge of the Galaxy for this line-of-sight. The line-of-sight, however, borders the giant HII region 
W5 in the Perseus Arm. We bracket the pulsar distance as being between W5 (2 kpc) and the 
Galaxy edge (9 kpc). 



PSR J0631+1036 - The DM = 125.3 cm" 3 pc (IZepka et al.lll996h is large for a source in the 
direction of the Galactic anticenter. The dark cloud LDN 1605, part of the active star-forming 
region 3 Mon, is in the line-of-sight. The pulsar could be inside the cloud, at ~0.75 kpc. Ionized 
material in the cloud could cause NE2001 to overestimate the distance. 

PSR Jl 124 — 5916 -It lies toward s the Carina arm where NE2001 biases are acute. The DM dis- 
tance is 5.7 kpc (jCamilo et al.ll2002a ). The kinematic distance of the associated SNR (G292.0+1.8) 
in dicates a lower limit of 6.2zfc0 .9 kpc (IGaensler fc Wallace! 120031 ) . The value in Table [5] is derived 
by iGonzalez &; Safi-Harbl (|2003l ) linking the X-ray absorption column with the extinction along the 
pulsar direction. 



PSR J1418-6058 - This pulsar is likely associate d with the PWN G313 .3+0.1, near the 
Kookaburra complex. A nearby HII region is at 13.4 kpc (jCaswell &; Hayneslll987l ) but could easily 
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be in the background. Such a large distance implies an unreasonably larg e gamma-ray efficiency- 
Table [5] lists a crude estimate of the distance range with the lower limit (jYadigaroglu &: Romani 
19971 ) taking the pulsar to be rela ted to one of t he near objects (Clust 3, CI Lunga 2 or SNR 
G312.4— 0.4), a nd th e higher limit (|Ng et al.ll2005l ) determined by applying the relation found by 



Possenti et al 



bv lGotthelj (j200J). 



(2002) and the correlation between pulsar X-ray photon index and luminosity given 



PSR J1709- 4429 - The NE2001 DM distance is 2.3±0.7 kpc (jTavlor et alJll993h. Kinematic 



dista nces give upper and lower limits of 3.2±0.4 kpc and 2.4±0.6 k pc, respectively (IKoribalski et al 



19951 ). T he X-ray flux from the neutron star detected by Chandra (jRomani et al,ll2005l ) and XMM- 
Newton (jMcGowan et al.ll2004l ) is compatible with a distance of 1.4-2.0 kpc. We assume the range 
1.4-3.6 kpc. 

PSR J1747-2958 - The pulsar is as sociated with the P WN G359. 23-0.82. HI measurements 
yield a dista nce upper limit of 5 .5 kpc (jUchida et al.lll992l ). but the DM (101 pccm -3 ) suggests 
2.0±0.6kpc (|Camilo et al.ll2002bl ). The X-ray absorbing column detected by Chandra is between 4 
and 5 kpc, while the close r value of 2 kpc would imply that an otherwise unknown molecular cloud 
lies in front of the pulsar (jGaensler et al.ll2004l ). A range of 2-5 kpc is used in our analysis. 



P SR J2021+3651 - The DM distance of ~12 kpc implies a high gamma-ray conversion effi- 



ciency ((Roberts et al.ll2002l : lAbdo et al.ll2009ml ). The open cluster Berkeley 87 near the line-of-sight 
could be responsible for an electron column density higher than modeled by NE2001. The distance 
i n Table [5] comes from a Chandra X-ray observation of the pulsar and its surrounding nebula 
(jVan Etten et al.ll2008l ). A similar range (1.3-4.1 kpc) was obtained for the X-ray flux detected 
from the associated PWN. 

PSR J2032+4127 - The DM value (115 pc cm" 3 ) gives an NE2001 distance of 3.6 kpc. If 
the pulsar belongs to the star cluster Cyg OB2, it would be located at approximately 1.6 kpc 
(jCamilo et al.ll2009i ). In this text we use a range of 1.6-3.6 kpc for this source. 



PSR J2229+6114 - The distance derived from the X-ray absorption is ~3 kpc (jHalpern et al 



2001bl ). between the values from the DM (6.5 kpc; lHalpern et al.ll2001al ) and from the kinematic 
method (0.8 kpc: kothes et alihoOll ). 



Figure [3] shows a polar view of the distribution of known pulsars over the Galactic plane. When 
two different distances are listed, we plot the closer one. 



3.2. Spatial Distributions, Luminosity, and Other Pulse Properties 

Figure Q] shows the pulsars projected on the sky. A Gaussian fit to the Galactic latitude 
distribution for those with |fe| < 10° and having distance estimates yields a standard deviation of 
<7b = 3.5 ± 0.8 degrees. The distances range from d = 0.25 to 5.6 kpc, and we can use the most 
distant to place an upper limit on the scale height, obtaining h < <isin3?5 = 340 pc, close to the 
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typical scale height for all radio pulsars (the value used in lLvne et al.l (|l998l ) is 450 pc) 



The light curve peak separations A and the radio lags S from Table [3] are summarized in Figure 
UJ As we will discuss in Section 5, outer magnetosphere emission models predict correlations between 
these parameters. Figure [5] shows -Blc versus the characteristic age (r c ). The magnetic fields at 
the light cylinder for the detected MSPs are comparable to those of the other gamma-ray pulsars, 
suggesting that the emission mechanism for the two families may be similar. 

Table [4] lists L 7 and r], while Figure [6] plots L 7 vs. E. The dashed line indicates L 7 = E, while 
the dot-dashed line indicates L 7 oc E 1 / 2 , where 



L 7 = 4ird 2 f n G 



100- 



(7) 



The flux correction factor (jWatters et al.ll2009l ) is model-dependent and depends on the magnetic 
inclination and observer angles a and (. Both the outer gap and slot gap mode ls predict fa ~ 1, 



in contrast to earlier use of fn = 1/47T ~ 0.08 (in e.g. iThompson et al.lll994l ). or = 0.5 for 
MSPs (in e.g. iFierro et al.lll995l ). For simplicity, we use /n = 1 throughout the paper, which 
presumably induces an artificial spread in the quoted L 7 values. However, it is the quadratic 
distance dependence for L 7 that dominates the uncertainty in L 7 in nearly all cases. 

Gamma-rays dominate the total power L tot radiated by most known high-energy pulsars, that 
is, L to t ~ L~. The Crab is a notable exception, with X-ray luminosity Lx ~ 10L 7 . In Figure [6] we 



plot both L 7 and Lx + L^. Lx for E < 100 MeV is taken from Figure 9 of iKuiper et al.l (|200ll ). 



3.3. Associations 



Table [6] provides some alternate names and positional associations of the pulsars in this catalog 
with other astrophysical sources. For the EGRET 3EG, EGR, and GEV and AGILE AGL catalogs, 
the uncertainties in the localization of the counterparts is worse than for the LAT sources. In these 
cases, we consider a source is a possible counterpart to a LAT pulsar when the separation between 
the two positions is less than the quadratic sum of their 95% confidence error radii. 

We see that 25 of the 46 pulsars are associated with sources in the 3EG, EGR and GEV 
catalogs of EGRET sources, though 19 were seen only as unidentified unpulsed sources. A number 
of these uni dentified EGRET sources had previously b een associated with SNRs, PWNe, or other 
objects (e.g. IWalker et al.ll2003l : iDe Becker et al.ll2005l ). In all cases, the gamma-ray emission seen 
with the LAT is dominated by the pulsed emission. Of the 25 EGRET sources, 14 are gamma-ray- 
selected pulsars, an d 11 are radio-selected, including 2 MSPs. All 6 high-confidence EGRET pulsars 



( Nolan et al.lll996T) are dete c ted, and the 3 marginal EGRET det ections are confirmed as pulsars 



Ramanamurthv et al.lll996l ; iKaspi et al.l 12000 ; IKuiper et al.ll2000l ). The 21 sources without 3EG 



EGR, or GEV counterparts include 18 previously detected radio pulsars (6 of which are MSPs) and 
3 gamma-ray selected pulsars. 
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Not surprisingly, many of the young pulsars have S NR or PWN associations. At le ast 19 of 
the 46 pulsars are associated with a PWN and/or SNR ([Roberts et alj|2005l : iGreenl |2009i ) . We do 
not test here whether the gamma-ray flux from any of these pulsars includes a non-magnetospheric 
component, as might be indicated by spatially extended emission or a spectrum at pulse minimum 
not characteristic of a pulsar. Such studies are underway. 

At least 12 of the pulsars are associated with TeV sources, 9 of which are also associated with 
PWNe. Those pulsars with both TeV and PWN associations are typically young, with ages less 
than 20 kyr. 



4. Pulsar Flux Sensitivity 

In order to interpret the population of gamma-ray pulsars discovered with the LAT, we need 
to evaluate the sensitivity of our searches for pulsed emission. While the precise sensitivity at any 
location is a function of the local background flux, the pulsar spectrum, and the pulse shape, we can 
derive an approximate pulsed sensitivity by calculating the unpulsed flux sensitivity for a typical 
pulsar spectrum at all locations in the sky and correlating with the observed Zf test statistic for 
the ensemble of detected pulsars. 

Figures [7] and show the distributions of the cutoff energy and the photon index, respectively, 
for all the LAT-detected pulsars. The distribution of photon indices peaks in the range T = 1 — 2, 
and the distribution of cutoffs peaks at -^cutoff = 1 — 3 GeV. For a typical spectrum, we used T = 1.4 
and -Ecutoff = 2.2 GeV, values approximately equal to their respective weighted averages. 

We then generated a sensitivity map for unpulsed emission for the six-month data set used 
here. For each (l,b) location in the sky, we computed the DC flux sensitivity at a threshold 
likelihood test statistic TS = 25 integrated above 100 MeV, assuming the typical pulsar spectrum 
within the source P SF and an underlying diffuse gamma-ray flux from the rings_Galaxy_vO model 



( Abdo et al.ll2009a ). This is an earlier version of the publicly available 2 model gll_iem_v02, similar 
to the model used for the spectral analysis. We note that the likelihood calculation assumes that the 
source flux is small compared to the diffuse background flux within the PSF, which is appropriate 
for a source just at the detection limit. Finally, we converted this map to pulsed sensitivity by a 
simple scale factor that accounts for the correspondence between the Z\ periodicity test confidence 
level and the unpulsed likelihood TS for the detected pulsars. 

The resulting 5<r sensitivity map for pulsed emission is shown in Figure [9j Comparing the 
measured fluxes with the predicted sensitivities at the pulsar locations (Figure [10]) , we see that this 
5c7 limit indeed provides a reasonable lower envelope to the pulsed detections in this catalog. Thus 
the effective sensitivity for high latitude (e.g. millisecond) pulsars with known rotation ephemerides 
is 1 — 2 x 10~ 8 cm" 2 s _1 ; at low latitude there is large variation, with typical detection thresholds 
3 — 5x higher. We expect the threshold to be somewhat higher for pulsars found in blind period 
searches. Figure [TUl suggests that this threshold is 2 — 3x higher than that for pulsars discovered 
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in folding searches, with resulting values as high as 2 x 10 7 cm 2 s 1 on the Galactic plane. 

The LogN-LogS plot is shown in Figure QTJ The approximate N oc 1/S dependence expected 
for a disk population is apparent for the higher flux objects. This shows that while radio-selected 
pulsars are detected down to a threshold of 2 x 10~ 8 cm" 2 s , the faintest gamma-ray-selected 
pulsar detected has a flux ~ 3x higher at 6 x 10~ 8 cm~ 2 s^ 1 . It is interesting to note that, aside 
from the lower flux threshold for the former, the radio-selected and gamma-ray-selected histograms 
are well matched, suggesting similar underlying populations. 



5. Discussion 



The striking results of the early Fermi pulsar discoveries demonstrate the LAT's excellent 
power for pulsed gamma-ray detection. By increasing the gamma-ray pulsar sample size by nearly 
an order of magnitude and by firmly establishing the gamma-ray-selected (radio-quiet Geminga- 
type) and millisecond gamma-ray pulsar populations, we have promoted GeV pulsar astronomy to 
a major probe of the energetic pulsar population and its magnetospheric physics. Our large pulsar 
sample allows us both to establish patterns in the pulse emission possibly pointing to a common 
origin of pulsar gamma-rays and to find anomalous systems that may point to exceptional pulsar 
geometries and/or unusual emission physics. In this Section we discuss some initial conclusions 
drawn from the sample, recognizing that the full exploitation of these new results will flow from 
the detailed population and emission physics studies to follow. 



5.1. Pulsar Detectability 

A w idely-cited predict or of gamma-ray pulsar det ectabil i ty is t he spin-do wn flux at Earth E/d 2 
(see e.g. ISmith et al.ll2008l ). However, as argued by lArond (120061 ) (see also lHarding &; Muslimov 
2002! ) . it is natural in many models for the gamma-ray emitting gap to m aintain a fixed voltage 
drop. This implies that L 7 is simply proportional to the particle current (|Hardingill98ll ). which 
gives L 7 oc E 1 / 2 , i.e. gamma-ray efficiency increases with decreasing spin-down power down to 
E ~ 10 34 — 10 35 erg s _1 where the gap saturates at large efficiency. In Figure [T2l we show how our 
detected pulsars rank in E 1 / 2 /d 2 against the set of searched pulsars. We see that for both MSPs 
and young pulsars, the detected objects have among the largest values of this metric. The presence 
of missing objects among the detected pulsars is interesting, but must be treated with caution, 
as the detectability metric may be inflated by poor DM distances, or the sensitivity of the pulse 
search might be anomalously low due to high local background or unfavorable pulsar spectrum 
or pulse profile. Alternatively, some missing objects may be truly gamma-ray faint for the Earth 
line-of-sight. A more complete study of the implications of the pulsar non-detections and upper 
limits is in progress. 



To study the luminosity evolution in the observed pulsar population, we plot in Figure [6] our 
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present best estimate of L 7 against E, based on the pulsed flux measured for each pulsar. Two 
important caveats must be emphasized here. First, the inferred luminosities are quadratically 
sensitive to the often large distance uncertainties. Indeed, for many radio selected pulsars (green 
points) we have only DM-based distance estimates. For many gamma-ray-selected pulsars we have 
only rather tenuous SNR or birth cluster associations with rough distance bounds. Only a handful 
of pulsars have secure parallax-based distances. Second, we have assumed here uniform phase- 
averaged beaming across the sky ( fo= l). This is not realized for many emission models, especially 



for low E pulsars (jWatters et al.l 12009 ) . 



To guide the eye, Figure [6] shows lines for 100% conversion efficiency (L 7 = E) and a heuristic 
constant voltage line L 7 = (10 33 erg s~ l E) 1 / 2 . In view of the large luminosity uncertainties, we 
must conclude that it is not yet possible to test the details of the luminosity evolution. However, 
some trends are apparent and individual objects highlight possible complicating factors. For the 
highest E pulsars, there does seem to be rough agreement with the E 1 / 2 trend. However, large 
variance between different distance estimates for the Vela-like PSRs J2021+3651 and J1709— 4429 
complicate the interpretation. In the range 10 35 erg s~ x < E < 10 36 5 erg s" 1 , the L 7 seems 
nearly constant, although the lack of precise distance measurements limits our ability to draw 
conclusions. For example, the very large nominal DM distance of PSR J0248+6021 would require 
> 100% efficiency, and so is unlikely to be correct. Two other pulsars with apparent high efficiency 
(J1836+5925 and J2021+4026) are plotted including relatively bright unpulsed emission; this may 
be magnetospheric, but may also be a surrounding or nearby source. The association distances 
for the gamma-ray-selected pulsars must additionally be treated with caution. For example PSR 
J2021+4026 has a r c ~ 10 x larger than the age of the putative associated SNR 7 Cygni. Improved 
distance estimates in this range are the key to probing luminosity evolution. 

From 10 34 erg s^ 1 < E < 10 35 erg s 1 we have several nearby pulsars with reasonably ac- 
curate parallax distance estimates. However we see a wide range of gamma-ray efficiencies. This 
is the range over which, for both slot gap and outer gap models, the gap is expected to 'satu- 
rate' and use most of the a vailable potential to maintain the pair cascade. In slot gap models 



( Muslimov fc Harding! I2003I ). the break occurs at about 10 35 erg s , when the gap is limited by 



screening of the accelerating fiel d by pairs. The e fficiency below this saturation is predicted to 



be ~ 10%. In outer gap models (jZhang et al.1 12004 ). the break is predicted to occur at somewhat 
lower E ~ 10 34 erg s _1 . With the present statistics and uncertainties, it is not possible to discrim- 
inate between these model predictions except to note that both are consistent with the observed 
results. In some models the gap saturation dramatically affects the shape of the beam o n the sky 



and accordingly the flux conversion factor /q; for outer gap models IWatters et al.l (|2009l ) estimate 
fn ~ 0.1 — 0.15 for Geminga (similar values are obtained for J1836+5925), driving down the rather 
high inferred luminosity of these pulsars by an order of magnitude. In contrast, another pulsar 
with an accurate parallax distance, PSR J0659+1414, has an inferred luminosity 30 x lower than 
the E 1 / 2 prediction. Clearly, some parameter in addition to E controls the observed L 7 . Finally, 
for < 10 erg s" 1 the sample is dominated by the MSPs. These nearby, low luminosity objects 



- 24 - 



clearly lie below the E 1 ! 2 trend, and in fact seem more consistent with L 7 oc E. 

Upper limits on radio pulsars with high values of the spin-down flux at Earth or large E 1 / 2 /d 2 
can help constrain viable efficiency models. In practice, the modest present exposure, the large 
background in the Galactic plane and the need to rely on uncertain dispersion-based distance 
estimates limit the value of such constraints. Still, a few pulsars are already interesting; for example, 
using the DM-based distance, the sensitivity in Figure[9l and an assumed fa = 1, we find that PSR 
J1740+1000 shows less than 1/5 of the flux expected from the E 1 ! 2 (constant voltage) line in Figure 
[6l Similarly, PSRs J1357— 6429 and J1930+1852 have upper limits just below the expected fluxes. 
Further, some detected pulsars, e.g. PSRs J0659+1414 and J0205+6449, lie significantly below 
the constant voltage trend. We expect that as LAT exposure and the significance of such limits 
increase, we should obtain additional constraints on the factors controlling pulsar detectability. 

One likely can didate for the additio nal factor affecting ga mma-ray detectabilit y is beaming. For 
PSR J1930+1852 (jCamilo et alJl2002ah . X-ray torus fitting \Ns & Romanilbood ) suggests a small 
viewing angle |C| ~ 33°. In outer gap models this makes it highly unlikely that the pulsar will pro- 
duce strong emission on the Earth line-of-sight. Similarly it h as been argued that PSR J0659+1414 
has a small viewing angle £ < 20° ([Everett fc Weisberal200ll ) (but see IWeltevrede Wrightl (j2009l ) 
for a discussion of uncertainties). Again, strong emission from above the null charge surface is not 
expected for this £. One possible interpretation is that we are seeing slot gap or even polar cap 
emission from this pulsar, which is expected at this £. The unusual pulse profile and spectrum of 
this pulsar may allow us to test this idea of alternate emission zones. 

In discussing non-detections, we should also note that the only binary pulsar systems reported 
in this paper are the radio-timed MSPs. In particular, our blind searches are not, as yet, sensitive 
to pulsars that are undergoing strong acceleration in bin ary systems. However, we do expect such 
objects to exist. Population syntheses (jPfahl et al.ll2002l ) suggest that several percent of the young 
pulsars are born while retained in massive star binary systems. A few such systems are known in 
the radio pulsar sample (e.g. the TeV-detected PSR B1259— 63); we expect that with the gamma- 
ray signal immune to dispersion effects an appreciable number of pulsar massive-star binaries will 
eventually be discovered. In deed, it is enti rely possible that the bright gamma-ray binaries LSI 
+61° 303 (|Abdo et al.ll2009hh and LS 5039 dAbdo et al.ll2009ih may host pulsed GeV signals that 
have not yet been found. 



5.2. Pulsar Population 



With the above caveats about missing binary systems in mind, we can already draw some 
conclusions about the single gamma-ray pulsar population. For example, there are 17 gamma-ray 
selected pulsars with a faintest flux of ~ 6 x 10 _8 cm _2 s _1 ; there are 16 non-millisecond radio- 
selected pulsars to t his flux limit. Of course, some gamma-ray-selected objects can indeed be 
detected in the radio (jCamilo et al.ll2009l ). Indeed, the detection of PSR J1741— 2054 at Li.4GHz 1=3 
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0.03 mJy kpc 2 underlines the fact that the radio emission can be very faint. Deep searches for 
additional radio counterparts are underway. However, wit h deep radio observat i ons of several 



objects, e.g . Geminga, PSR J0007+7303, PSR J1836+5925 (jKassim fc Laziol 11999 : lHalpern et al 



20041 . 120071 ). providing no convincing detections, it is clear that some objects are truly radio faint. 
The substantial number of radio faint objects suggests that gamma-ray emission has an appreciably 
larger extent than the radio beams, such as expected in the outer gap (OG) and slot-gap/two pole 
caustic (SG/TPC) models. 

Population synthesis studies for normal (non-millisecond) pulsars predicted that LAT would 
detec t from 40-80 radio loud pulsars and com parable numbers of radio quiet pulsars in the first 
year (jGonthier et al,l [2004J; IZhang et al.l 120071 ) . The ratio of radio-selected to gamma-ray-selected 
gamma-ray pulsars has been noted as a particularly sensitive discriminator of mo dels, since the 



outer magnetosphere models predict much smaller rati os than polar cap models (jHarding et al 



20071 ). Studies of the MSP population (jStorv et al.ll2007l ) predicted that LAT would detect around 
12 radio-selected and 33-40 gamma-ray-selected MSPs in the first year, in rough agreement with the 
number of radio-selected MSPs seen to date (searches for gamma-ray selected MSPs have not yet 
been conducted). Thus, in the first six months the numbers of LAT pulsar detections are consistent 
with the predicted range, and the large number of gamma-ray selected pulsars discovered so early 
in the mission points towards the outer magnetosphere models. 

We can in fact use our sample of detected gamma-ray pulsars to estimate the Galactic birthrates. 
For each object with an available distance estimate, we compute the maximum distance for detec- 
tion from D max = Z?est(-^7/-^min) 1 ^ 2 ) where D est comes from Table [SJ the photon flux Fioo from 
Table H and F min from Figure M We limit D 

max to 15 kpc, and compare V , the volume enclosed 
within the estimated source distance, to V m3i ^, that enclosed within the maximum distance, for a 
Galactic disk with radius 10 kpc and thickness lkpc. If we assume a blind search threshold 2x 
higher than that for a folding search at a given sky position, the inferred values of (V/V max ) are 
0.49, 0.59 and 0.55 for the radio-selected young pulsars, millisecond pulsars and gamma-ray-selected 
pulsars, respe ctively. These are close to 0.5, the value expected for a population uniformly filling a 
given volume (|Schmidtlll968l ): the MSP value is somewhat high as our sample includes three objects 
detected at < 5a. The value for the gamma-ray-selected pulsars is also high but is controlled by 
PSR J2032+4127. If we exclude this object from the sample, we get {V/V max ) = 0.5 at an effective 
threshold of 3x the ephemeris- folding value. 

Although we do not attempt a full population synthesis here, if we assume that the pulsar 
characteristic ages are the true ages, our sample can give rough estimates for local volume birthrates: 
8 x 10 kpc" 3 yr" 1 (young radio-selected), 4 x 10 -5 kpc~ 3 yr _1 (young gamma-ray-selected, 2x 
threshold) and 2 x 10 -8 kpc" 3 yr" 1 (MSP). Note that only half of the gamma-ray-selected objects 
have distance estimates. If we assume that the set without distance information has comparable 
luminosity, the gamma-ray-selected birthrate is ~ 2x larger. These estimates retain appreciable 
uncertainty; for example if the effective blind search detection threshold is 3 x that for folding, the 
inferred gamma-ray-selected birthrate increases by an additional ~ 65%. If we extrapolate these 
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local birthrates to a full disk with an effective radius of lOkpc we get 1/100 yr (radio-selected 
young pulsars), 1/100 yr (gamma-ray-selected pulsars) and 1/(6 x 10 5 yr) (radio-selected MSPs). 
Normally in estimating radio pulsar birthrates one would correct for the radio beaming fraction. 
However if young gamma-ray-selected pulsars are simply similar objects viewed outside of the radio 
beam, this would result in double-counting. In any case one infers a total Galactic birthrate for 
energetic pulsars of ~ 1/50 yr, with gamma-ray-selected objects representing half. This represents a 
large fraction of the estimated Galactic supernova rate, so clearly more careful population syntheses 
will be needed to see if these numbers are compatible. 



5.3. Trends in Light Curves and Other Observables 

The pulse shape properties can also help us to probe the geometry and physics of the emission 
region. The great majority of the pulsars show two dominant, relatively sharp peaks, suggesting 
that we are seeing caustics from the edge of a hollow cone. When a single peak is seen, it tends to 
be broader, suggesting a tangential cut through an emission cone. This picture is realized in the 
OG and the high altitude portion of the SG models. 

For the radio-emitting pulsars, we can compare the phase lag between the radio and first 
gamma-ray peak 5 wit h the separation of the two gamma-ray peaks A. As first pointed out in 



Romani fc Yadigaroglul (|1995l ). these should be correlated in outer magnetosphere models — this 



is indeed seen (Figu re HI). The dis t ribut ion can be compared with predictions of the TPC and 



OG models shown in lWatters et al.l (|2009l ). The 5 — A distribution and in particular the presence 
of A ~ 0.2 — 0.3 values appear to favor the OG picture. However, there are a greater number 
having A ~ 0.4 — 0.5, which favors TPC models. In Figure [TBI we show the peak separation as a 
function of pulsar spin-down luminosity — the A distribution appears to be bimodal, with no strong 
dependence on pulsar E (or age). A full comparison will require detailed population models, which 
are being created. It may also be hoped that the precise distribution of measured values can help 
probe details of the emission geometry. In particular, whenever we have external constraints on the 
viewing angle £ (typically from X-ray images of the PWN) or magnetic inclination a (occasionally 
measured from radio polarization), then the observed values of 5 and A become a powerful probe 
of the precise location of the emission sheet within the magnetosphere. This can be sensitive to the 
field perturbations from magnetospheric currents and hence can probe the global electrodynamics 
of the pulsar magnetosphere. 

If one examines the energy dependence of the light curves of both the radio-selected and 
gamma-ray-selected pulsars, a decrease in the P1/P2 ratio with increasing energy seems to be 
a common feature. However, the P1/P2 ratio evolution does not occur for all pulsars, notably 
J0633+0632, J1028-5819, J1124-5916, J1813-1246, J1826-1256, J1836+5925, J2021+3651, and 
J2238+59. Most of these pulsars have two peaks with phase separation of ~0.5 and little or no 
bridge emission between the peaks. Perhaps the lack of P1/P2 energy evolution is connected with 
an overall symmetry of the light curve. 
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The LAT pulsar sample also shows evidence of trends in other observables that may offer 
additional clues to the pulsar physics. While the detected objects have a wide range of surface 
magnetic fields, their inferred light cylinder magnetic fields B^c are uniformly relatively large 
(> 10 3 G). Indeed, the LAT-detected MSPs are those with the highest light cylinder fields with 
values very similar to those of the detected normal pulsars. Comparison of the spectral cut-off 
^cutoff with surface magnetic field shows no significant correlation. This evidence argues against 
classical low altitude polar cap models supported by 7-B cascades. However, there is a weak 
correlation of Cutoff with £?lc > as shown in Figure d It is interesting that the values of -E C utoff 
have a range of only about a decade, from 1 to 10 GeV, and that all the different types of pulsars 
seem to follow the same correlation. This strongly implies that the gamma-ray emission originates 
in similar locations in the magnetosphere relative to the light cylinder. Such a correlation of .E cu toff 
with i?LC is actually expected in all outer magnetosphere models where the gamma-ray emission 
primarily comes from curvature radiation of electrons whose acceleration is balanced by radiation 
losses. In this case, 

'£„\ 3 / 4 
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in m c c 2 , where A c is the electron Compton wavelength, E\\ is the electric field that accelerates 
part icles parallel to the magne tic field a nd p c - is the magnetic field radius of curvature. In both 



SG dMuslimov & Hardindl2004h and OG (jZhang et al.1 12004 : lHirotanill2008h models, E\\ oc B LC w 



where w is the gap width. All these models give values of E cuto g that are roughly consistent with 
those measured for the LAT pulsars. Although p c ~ Rlc, the gap widths are expected to decrease 
with increasing B^c > so that E cuto g is predicted to be only weakly dependent on £>lc in most outer 
magnetosphere models, as observed. 

The detection of pulsed flux at E max > a few GeV provides additional, physical motivation for 
high altitude emission, since one expects strong (hyper-exponential) attenuation from 7-B — > e + e~ 
absorption in the high magnetic fields at low altitudes of a near-surface po lar gap. For a polar 
cap model with spin period P and surface field 10 12 i?i2G, Equation 1 of iBarina ( 2004 ) gives 
r ^ (-E'max-Bi2/l-76GeV) 2 / 7 P" 1 / 7 10 6 cm. While the largest minimum r is derived from the 25 GeV 
detection of the Crab by MAGIC ([Albert et al.ll2008l ) , significant minimum altitudes of 2 to 3 stellar 
radii are found for many LAT pulsars with large Byi and high E mto ff, assuming maximum energies 
set at around 2.5-E 0u t o ff ■ Such altitudes are inconsistent with the r < 1 — 2 stellar radii of polar cap 
models, hence implicating outer magnetosphere (e.g. TPC or OG) models where the bulk of the 
emission occurs at tens to hundreds of stellar radii. 

In Figure [HJ we see a general trend for the young pulsars to show a softer spectrum a t larg e 
E, although there is a great deal of scatter; a similar trend was noted in iThompson et al.l ()19 99). 
This may be indicative of higher pair multiplicity, which would steepen the spectrum for the more 
energetic pulsars, either by steepening the spectrum of the curvature radiation-generating primary 
electrons (|Romanilll996l) or by inclusion of an additional soft sp ectral component associated with 
robust pair formation ( Takata Sz Changll2007 : Harding et al.12008 ). In either case, one would expect 
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steepening from the simple monoenergetic curvature radiation spectrum T = 2/3 for the higher E 
pulsars. Interestingly, the MSPs do not extend the trend to lower E. Of course EGRET (and 
now the LAT) find strong variations of photon index with phase for the brighter pulsars. A full 
understanding of photon index trends will doubtless require phase-resolved modeling. 



6. Conclusion 

The new gamma-ray pulsar populations established by early LAT observations show that we 
are detecting many nearby young pulsars. In addition we are detecting the millisecond pulsars 
with the highest spin-down flux at Earth. Thus we see that the LAT is providing a new, local, 
but relatively unbiased view of the energetic pulsar population (see Figures [H [2j and [3]). These 
detections provide a new window into pulsar demographics and physics. 

We conclude that a large fraction of the local energetic pulsars are GeV emitters. There is also 
a significant correlation with X-ray and TeV bright pulsar wind nebulae. Conversely, we have now 
uncovered the pulsar origin of a large fraction of the bright unidentified Galactic EGRE T sources, 



as proposed by several authors (jKaaret fe Cottaml Il996l : lYadigaroglu k, Romanil 1 1 99 71 ) . We have 



also found plausible pulsar counterparts for several previously detected TeV sources. In this sense 
the "mystery" of the unidentified EGRET sources is largely solved. It is possible that the two 
LAT-detected massive binaries (LSI +61° 303, LS 5059) and some of the remaining unidentified 
sources also contain spin-powered pulsars. Thus we expect that the LAT pulsar population will 
increase, with both the detection of binary gamma-ray pulsars and fainter and more distant pulsars. 

The light curve and spectral evidence summarized above suggests that these pulsars have high 
altitude emission zones whose fan-like beams scan over a large portion of the celestial sphere. This 
means that they should provide a relatively unbiased census of energetic neutron star formation. 
A rough estimate of the young gamma-ray pulsar birthrate extrapolating from our local sample 
suggests a Galactic birthrate as high as ~ 1/50 yr, a large fraction of the estimated Galactic 
supernova rate. Gamma-ray detectable MSPs in the Galactic field are born rarely, ~ 1/6 x 10 5 yr, 
but with their long lifetimes are inferred to contribute comparably to the number of (in principle) 
detectable Galactic gamma-ray pulsars. 

The data also advance our understanding of emission zone physics. It is now clear that the 
gamma-ray emission from the brightest pulsars arises largely in the outer magnetosphere. The pho- 
ton emission also accounts for a large fraction of the spin-down luminosity, increasing as the pulsars 
approach E ~ io 33 ~ 34 erg s" 1 . While these wide, bright beams are a boon for population studies, as 
noted above, they represent a challenge for theorists trying to understand pulsar magnetospheres. 
Further LAT pulsar observations and, in particular, the high quality, highly phase-resolved spectra 
now being obtained for the brightest LAT pulsars will surely sharpen this challenge. 
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Table 1. Measured and intrinsic parameters of LAT-detected pulsars 
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r i 


201.1 


8.3 


385 


55.0 


110 


3.8 


0.7 


3.7 


J0742— 2822 


r ; 


243.8 


—2.4 


167 


16.8 


160 


14.3 


3.3 


15 


707514-1807 

'J U i Ul | J.UU 1 




202.7 


21.1 


3.5 


6xl0 -6 


8. Ox 10 6 


0.6 


32.3 


3.2 


J0835— 4510 




263.6 


—2.8 


89.3 


124 


11 


688 


43.4 


1100 


71028—5819 


r 1 


285.1 


—0.5 


91.4 


16.1 


90 


83.2 


14.6 


0.36 


J1048— 5832 


j, m 


287.4 


0.6 


124 


96.3 


20 


201 


16.8 


6.5 


J1057— 5226 


j. n 


286.0 


6.6 


197 


5.8 


540 


3.0 


1.3 


11 


J1124— 5916 


2' 


292.0 


1.8 


135 


747 


3 


1190 


37.3 


0.08 


J1418— 6058 


g 6 
& 


313.3 


0.1 


111 


170 


10 


495 


29.4 


2 ' 3 < 0.06 


J1420— 6048 


r * 


313.5 


0.2 


68.2 


83.2 


13 


1000 


69.1 


0.9 


J1459— 60 


„■ b 
& 


317.9 


— 1.8 


103 


25.5 


64 


91.9 


13.6 


2 < 0.2 


J1509— 5850 


r * 


320.0 


—0.6 


88.9 


9.2 


150 


51.5 


11.8 


0.15 


J1614— 2230 


mb d 


352.5 


20.3 


3.2 


4x 10~ 6 


1.2xl0 6 


0.5 


33.7 




J1709— 4429 




343.1 


—2.7 


102 


93.0 


18 


341 


26.4 


7.3 


J1718— 3825 


r i 


349.0 


—0.4 


74.7 


13.2 


90 


125 


21.9 


1.3 


J1732— 31 


s 6 

6 


356.2 


0.9 


197 


26.1 


120 


13.6 


2.7 


2 < 0.2 


71 741 —2054 


a- M 
6 


6.4 


4.6 


414 


16.9 


390 


0.9 


0.3 


2 0.16 


J1744-1134 


m d 


14.8 


9.2 


4.1 


7xl0~ e 


9xl0 6 


0.4 


24.0 


3 


J1747-2958 


r ° 


359.3 


-0.8 


98.8 


61.3 


26 


251 


23.5 


0.08 


J1809-2332 


g 6 


7.4 


-2.0 


147 


34.4 


68 


43.0 


6.5 


2 - 3 < 0.06 


J1813-1246 


g " 


17.2 


2.4 


48.1 


17.6 


13 


626 


76.2 


2 < 0.2 


J1826-1256 




18.5 


-0.4 


110 


121 


11 


358 


25.2 


2 - 3 < 0.06 


J1833-1034 


r ° 


21.5 


-0.9 


61.9 


202 


5 


3370 


137.3 


0.07 


J1836+5925 




88.9 


25.0 


173 


1.5 


1800 


1.2 


0.9 


4 < 0.007 


J19074-06 


g i>,r 


40.2 


-0.9 


107 


87.3 


19 


284 


23.2 


< 0.02 


J1952+3252 


r n 


68.8 


2.8 


39.5 


5.8 


110 


374 


71.6 


1 


J1958+2846 


g " 


65.9 


-0.2 


290 


222 


21 


35.8 


3.0 




J20214-3651 


r p 


75.2 


0.1 


104 


95.6 


17 


338 


26.0 


0.1 


J20214-4026 


g M 


78.2 


2.1 


265 


54.8 


77 


11.6 


1.9 




J20324-4127 


g M 


80.2 


1.0 


143 


19.6 


120 


26.3 


5.3 


2 0.24 


J2043+2740 


r « 


70.6 


-9.2 


96.1 


1.3 


1200 


5.6 


3.6 


5 3 


J2124-3358 


m d 


10.9 


-45.4 


4.9 


12xl0~ 6 


0.6xl0 6 


0.4 


18.8 


1.6 


J22294-6114 


r m 


106.6 


2.9 


51.6 


78.3 


11 


2250 


134.5 


0.25 
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Table 1 — Continued 



PSR 


Type, 


1 


b 


P P age r c 


E B LC 


Suoo 




Ref. 


(°) 


(°) 


(ms) (10~ 15 ) (kyr) 


(10 34 erg s- 1 ) (kG) 


(mjy) 



J2238+59 g b 106.5 0.5 163 98.6 26 90.3 8.6 



Note. — The first two columns are pulsar names and types: r for radio-selected, g for gamma-ray- 
selected, m for MSPs, and b for binary pulsars. The 3 rd and 4 th columns are Galactic coordinates for 
each pulsar. The 5 th and 6 th columns list the period (P) and its first derivative (P), corrected for the 
Shklovskii effect (see text). Following are the characteristic age t c (column 7), the spin-down luminosity 
E (column 8), and the magnetic field at the light cylinder Blc (column 9) . The last column is the 
radio flux densit y at 1400 MHz, or an up per limit when one is available. These values are taken from the 
ATNF database ^Manchester et al.ll2005l) exce pt for the noted entries wh ere: ( l) jHalpern et aljlioplh ; (2) 
iCamilo et afll2009h ; (3) ^Roberts et aljbooi) ; (4) iHalpern et alj|2007l) ; (5) jptav et al.lll996Tl . Note that 
PSR J1509-5850 should not be confused with PSR B1509-58 observed by CGRO. 



I Cognard et al 



References. — References to Fermi L AT pub licatio ns specific t o these pulsars: a jAbdo et al.l 120081 1 
b jAbdo et al l |2009ch ; c jAbdo et al.l l2009nTl ; d jAbdo et al.l l2009tj) ; e jAbdo et al.l l2009otl : f 



'l201Ch ; h jAbdo et al.l l2010dl) ~ i jWeltevrede et al 



l Abdo et al.l2009e) : I jAbdo et al|2009eh : mjAbdo et al.l2009jl~n" Abdo et alj2010a) : o jCamilo et al.l 



120101) : v llAbdo et al.ll2009mh : iNoutsos et alj|2010h 
iCamilo et alibOOSh 



201Ch : 7 jAbdo et al.l l2009all 



jAbdo et al.ll2010bl' l : s llAbdo et al.lbOlOch 
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Table 2. Pulsation detection significances for LAT-detected pulsars 



PSR Z\ value H value maxROI(°) ObsID 



J0007+7303 


2072.1 


2371.8 


1.0 


L 


J0030+0451 


121.1 


362.7 


1.0 


N 


J0205+6449 


90.9 


206.0 


1.0 


G, J 


J0218+4232 


24.7 


22.5 


1.0 


N, W 


J0248+6021 


57.5 


75.1 


0.5 


N 


J0357+32 


422.7 


450.7 


1.0 


L 


J0437-4715 


126.9 


153.6 


1.0 


1' 


J0534+2200 


4397.8 


15285.0 


1.0 


N, J 


J0613-0200 


93.6 


139.9 


1.0 


N 


J0631+1036 


48.6 


44.8 


1.0 


N, J 


J0633+0632 


230.2 


573.3 


1.0 


L 


J0633+1746 


10053.6 


20346.4 


1.0 


L 


J0659+1414 


80.5 


99.0 


1.0 


N, J 


J0742-2822 


38.9 


44.9 


1.0 


N, J 


J0751+1807 


29.7 


26.5 


1.0 


N 


J0835-4510 


26903.9 


74716.7 


1.0 


P 


J1028-5819 


291.5 


915.9 


0.5 


1' 


J1048-5832 


208.5 


634.0 


1.0 


P 


J1057-5226 


1668.9 


1772.4 


1.0 


P 


J1124-5916 


93.5 


179.9 


1.0 


L 


J1418-6058 


230.1 


343.7 


1.0 


L 


J1420-6048 


104.7 


114.4 


1.0 


P 


J1459-60 


148.2 


159.3 


1.0 


L 


J1509-5850 


71.6 


73.3 


0.5 


1' 


J1614-2230 


36.2 


69.5 


0.5 


G 


J1709-4429 


4680.1 


5612.1 


1.0 


1' 


J1718-3825 


111.9 


109.8 


0.5 


N, P 


J1732-31 


141.2 


279.6 


1.0 


L 


J1741-2054 


332.6 


355.9 


1.0 


L 


J1744-1134 


28.4 


38.1 


1.0 


N 


J1747-2958 


47.2 


69.0 


0.5 


G 


J1809-2332 


589.3 


1562.5 


1.0 


L 


J1813-1246 


140.0 


162.0 


1.0 


L 


J1826-1256 


442.4 


979.0 


1.0 


L 


J1833-1034 


35.2 


87.6 


1.0 


G 


J1836+5925 


349.2 


385.3 


1.0 


L 


J1907+06 


257.1 


521.0 


1.0 


L 


J1952+3252 


464.8 


1008.8 


1.0 


J, N 


J1958+2846 


146.9 


233.1 


1.0 


L 


J2021+3651 


1433.5 


4603.7 


1.0 


G, A 


J2021+4026 


222.0 


275.8 


1.0 


L 


J2032+4127 


224.9 


485.9 


0.5 


L 


J2043+2740 


28.2 


38.2 


1.0 


N, J 


J2124-3358 


77.8 


80.9 


1.0 


N 


J2229+6114 


1026.0 


1237.4 


1.0 


G, J 


J2238+59 


135.8 


373.0 


1.0 


L 
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Table 2 — Continued 



PSR Zf value H value maxROI(°) ObsID 



No te. — C olumns 2 and 3 list the Z| iBuccheri et al.l 
ll983T) and H ide Jager et al.lll989h periodicity test values 
for E > 0.3 GeV respectively. Detection of gamma-ray 
pulsations are claimed when the significance of the period- 
icity test exceeds 5c, with the exceptions of J0218+4232, 
J0751+1807, J1744-1134, and J2043+2740 as described 
in Section 2.1.3. A significance greater than 5cr corre- 
sponds to Z| > 36 and H > 42 ; greater than 7cr corre- 
sponds to Z| > 61 ; and greater than 10a corresponds to 
Z| > 114 (see Section 2.1.3). Column 4 gives the maxi- 
mum angular radius (maxROI) around the pulsar position 
within which gamma-ray events were searched for pulsa- 
tions. The final column indicates the observatories that 
provided ephemerides (see Section 2.1.1 for details): "A" 
— Arecibo telescope; "G" - Green Bank Telescope; "J" - 
Lovell telescope at Jodrcll Bank; "L" - Large Area Tele- 
scope; "N" - Nancay Radio Telescope; "P" - Parkes radio 
telescope; "W" - Westerbork Synthesis Radio Telescope. 
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Table 3. Pulse shape parameters of LAT-detected pulsars 



PSR 


Type" 


Peak 


Radio lag 


7-ray peak separation Off-pulse definition 






multiplicity 


S 


A <j> 



J0007+7303 


g 


2 




0.23 ± 0.01 


0.29 


- 0.87 


J0030+0451 


m 


2 


0.18 ± 0.01 


0.44 ± 0.01 


0.68 


- 0.12 


J0205+6449 


r 


2 


0.08 ± 0.01 


0.50 ± 0.01 


0.64 


- 0.02 


J0218+4232 


111 


2 


0.32 ± 0.02 


0.36 ± 0.02 


0.84 


- 0.16 


J0248+6021 


r 


1 


0.35 ± 0.01 




0.71 


- 0.19 


J0357+32 


g 


1 






0.34 


- 0.86 


J0437-4715 


in 


1 


0.43 ± 0.02 




0.60 


- 0.20 


J0534+2200 


r 


2 


0.09 ± 0.01 


0.40 ± 0.01 


0.62 


- 0.98 


J0613-0200 


in 


1 


0.42 ± 0.01 




0.56 


- 0.16 


J0631+1036 


r 


1 


0.54 ± 0.02 




0.80 


- 0.20 


J0633+0632 


g 


2 




0.48 ± 0.01 


0.09 


- 0.45 


J0633+1746 


g 


2 




0.50 ± 0.01 


0.24 


- 0.54 


J0659+1414 


r 


1 


0.21 ± 0.01 




0.40 


- 1.00 


J0742-2822 


r 


1 


0.61 ± 0.02 




0.84 


- 0.44 


J0751+1807 


in 


1 


0.43 ± 0.02 




0.63 


- 0.99 


J0835-4510 


r 


2 


0.13 ± 0.01 


0.43 ± 0.01 


0.66 


- 0.06 


J1028-5819 


r 


2 


0.19 ± 0.01 


0.47 ± 0.01 


0.76 


- 0.12 


J1048-5832 


r 


2 


0.15 ± 0.01 


0.42 ± 0.02 


0.64 


- 0.04 


J1057-5226 


r 


2 


0.35 ± 0.05 


0.20 ± 0.07 


0.72 


- 0.20 


J1124-5916 


r 


2 


0.11 ± 0.01 


0.49 ± 0.01 


0.70 


- 0.06 


J1418-6058 


g 


2 




0.47 ± 0.01 


0.54 


- 0.90 


J1420-6048 


r 


2'' 


0.26 ± 0.02 


0.18 ± 0.02 


0.60 


- 0.10 


J1459-60 


g 


2 




0.15 ± 0.03 


0.34 


- 0.78 


J1509-5850 


r 


2'' 


0.18 ± 0.03 


0.20 ± 0.03 


0.52 


- 1.00 


J1614-2230 


in 


2 


0.19 ± 0.01 


0.51 ± 0.01 


0.92 


- 0.14 


J1709-4429 


r 


2 


0.24 ± 0.01 


0.25 ± 0.01 


0.66 


- 0.14 


J1718-3825 


r 


1 


0.42 ± 0.02 




0.68 


- 0.20 


J1732-31 


g 


2 




0.42 ± 0.02 


0.49 


- 0.93 


J1741-2054 


g 


2 


0.30 ± 0.01 


0.18 ± 0.02 


0.67 


- 0.19 


J1744-1134 


in 


1 


0.83 ± 0.02 




0.08 


- 0.44 


J1747-2958 


r 


2 


0.18 ± 0.01 


0.42 ± 0.04 


0.64 


- 0.10 


J1809-2332 


g 


2 




0.35 ± 0.01 


0.41 


- 0.89 


J1813-1246 


g 


2 




0.47 ± 0.02 


0.56 


- 0.90 


J1826-1256 


g 


2 




0.47 ± 0.01 


0.54 


- 0.94 


J1833-1034 


r 


2 


0.15 ± 0.01 


0.44 ± 0.01 


0.68 


- 0.10 


J1836+5925 


g 


2 




0.48 ± 0.01 






J1907+06 


g 


2 




0.40 ± 0.01 


0.46 


- 0.94 


J1952+3252 


r 


2 


0.15 ± 0.01 


0.49 ± 0.01 


0.68 


- 0.08 


J1958+2846 


g 


2 




0.45 ± 0.01 


0.55 


- 0.95 


J2021+3651 


r 


2 


0.17 ± 0.01 


0.47 ± 0.01 


0.70 


- 0.04 


J2021+4026 


g 


2 




0.48 ± 0.01 






J2032+4127 


g 


2 


0.15 ± 0.01 


0.50 ± 0.01 


0.60 


- 0.92 


J2043+2740 


r 


2 


0.20 ± 0.01 


0.36 ± 0.01 


0.64 


- 0.08 


J2124-3358 


in 


1 


0.86 ± 0.02 




0.92 


- 0.58 


J2229+6114 


r 


1 


0.49 ± 0.01 




0.64 


- 0.14 
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Table 3 — Continued 



PSR 


Type" 


Peak 
multiplicity 


Radio lag 
8 


7-ray peak separation 
A 


Off-pulse definition 


J2238+59 


g 


2 




0.50 ± 0.01 


0.60 - 0.92 



a Types are r=radio-selected, g=gamma-ray-selected, m=millisecond 

b For some pulse profiles the current dataset does not all ow clear discrimination b etween a single, 
broad pulse and two unresolved pulses. See the discussion in lWeltevrede et al. I tepid) regarding PSRs 
J1420-6048 and J1509-5850. 



Note. — Light curve shape parameters evaluated from the full energy range light curve (see Section 
2.1.3). These include the peak multiplicity (3rd column), the lag <5 of the first gamma peak from the 
main radio peak for the radio-detected pulsars (4th column), and the phase difference A between the 
main gamma-ray peaks (5th column). Column 6 lists the off-pulse phase range used in the spectral 
analysis. The boldface entries for PSR J 1124-5916 are the corrected values as per an Erratum sent 
to the ApJ in December 2010. 



Table 4. Spectral fitting results for LAT-detected pulsars 



PSR 


Type" 


Photon Flux (F 10 o) 
(lCT 8 ph cm- 2 s-') 


Energy Flux (Gioo) 
(10 -11 erg cm -2 s —1 ) 




r 




-^cutoff 

(GeV) 


TS 


TS cuto ff 


Luminosity 
(10 33 ergs" 1 ) 


Efficiency'' 
(/« = 1) 


J0007+7303 


g 


30.7 ± 1.3 


38.2 ± 1.3 


1 


.38 ± 


0.05 


4.6 ± 0.4 


7384 


274.7 


89 ± 38 


0.20 ± 0.08 


J0030+0451 


m 


5.83 ± 0.78 


5.26 ± 0.42 


1 


.22 ± 


0.19 


1.8 ± 0.4 


960 


59.2 


0.57 ± 0.35 


0.17 ± 0.10 


J0205+6449 


r 


13.2 ± 2.0 


6.64 ± 0.65 


2 


.09 ± 


0.17 


3.5 ± 1.4 


346 


12.5 


54 - 81 


0.002 - 0.003 


J0218+4232 


m 


6.2 ± 1.7 


3.62 ± 0.64 


2 


.02 ± 


0.28 


5.1 ± 4.2 


119 


4.7 


27 - 69 


0.11 - 0.29 


J0248+6021 


r 


3.7 ± 1.8 


3.07 ± 0.70 


1 


.15 ± 


0.59 


1.4 ± 0.6 


103 


18.5 


15 - 300 


0.07 - 1.4 


J0357+32 


g 


10.4 ± 1.2 


6.38 ± 0.44 


1 


.29 ± 


0.22 


0.9 ± 0.2 


949 


71.6 






J0437-4715 


m 


3.65 ± 0.84 


1.86 ± 0.26 


1 


.74 ± 


0.38 


1.3 ± 0.7 


172 


9.9 


0.054 ± 0.008 


0.02 ± 0.003 


J0534+2200 c 


r 


209 ± 4 


130.6 ± 3.4 


1 


.97 ± 


0.02 


5.8 ± 0.5 


21507 


80.2 


620 ± 310 


0.001 ± 0.001 


J0613-0200 
J0631+1036 


m 
r 


3.38 ± 0.85 
2.8 ± 1.2 


3.23 ± 0.42 
3.04 ± 0.61 


1 
1 


.38 ± 
.38 ± 


0.29 
0.42 


2.7 ± 1.0 

3.6 ± 1.8 


285 
86 


18.5 
10.0 


o.89 +82 

2.0 - 48 


07 + 06 
u - u( -0.03 

0.01 - 0.27 


J0633+0632 


g 


8.4 ± 1.4 


8.0 ± 0.77 


1 


.29 ± 


0.22 


2.2 ± 0.6 


370 


50.8 






J0633+1746 


g 


305.3 ± 3.5 


338.1 ± 3.5 


1 


.08 ± 


0.02 


1.9 ± 0.05 


62307 


5120.4 


25+ 24 


0-78 ±g;S 


J0659+1414 


r 


10 ± 1.4 


3.17 ± 0.36 


2 


.37 ± 


0.50 


0.7 ± 0.5 


206 


6.9 


0.31 ± 0.08 


0.01 ± 0.002 


J0742-2822 


r 


3.18 ± 1.2 


1.82 ± 0.42 


1 


.76 ± 


0.48 


2.0 ± 1.4 


47 


1.2 


9.0 ±l 2 


0-07 +°;g? 


J0751+1807 


m 


1.35 ± 0.66 


1.09 ± 0.38 


1 


.56 ± 


0.70 


3.0 ± 4.3 


37 


3.8 


0-47 +1 35 


0-08 ±g;SS 


J0835-4510 


r 


1061 ± 7.0 


879.4 ± 5.4 


1 


.57 ± 


0.01 


3.2 ± 0.1 


219585 


5971.0 


87 ± 12 


0.01 ± 0.002 


J1028-5819 


r 


19.6 ± 3.1 


17.7 ± 1.4 


1 


.25 ± 


0.20 


1.9 ± 0.5 


620 


75.1 


120 ± 73 


0.14 ± 0.09 


J1048-5832 


r 


19.7 ± 3.0 


17.2 ± 1.3 


1 


.31 ± 


0.18 


2.0 ± 0.4 


881 


81.8 


150 ± 90 


0.08 ± 0.05 


J1057-5226 


r 


30.45 ± 1.7 


27.2 ± 0.98 


1 


.06 ± 


0.10 


1.3 ± 0.1 


4961 


366.3 


17 ± 9 


0.56 ± 0.31 


J1124-5916 


r 


5.2 ± 1.8 


3.79 ± 0.70 


1 


.43 ± 


0.40 


1.7 ± 0.7 


111 


16.7 


100 tit 


01 +0.003 
u - ul -0.004 


J1418-6058 


g 


27.7 ± 8.3 


23.5 ± 3.8 


1 


.32 ± 


0.24 


1.9 ± 0.4 


162 


54.1 


110 - 700 


0.02 - 0.14 


J1420-6048 


r 


24.2 ± 7.9 


15.8 ± 3.5 


1 


.73 ± 


0.24 


2.7 ± 1.0 


63 


21.4 


590 ± 380 


0.06 ± 0.04 


J1459-60 


g 


17.8 ± 3.4 


10.56 ± 1.2 


1 


.83 ± 


0.24 


2.7 ± 1.1 


337 


21.1 






J1509-5850 


r 


8.7 ± 1.4 


9.7 ± 1.2 


1 


.36 ± 


0.28 


3.5 ± 1.1 


262 


26.3 


78 ± 49 


0.15 ± 0.10 


J1614-2230 


m 


2.89 ± 1.2 


2.74 ± 0.50 


1 


.34 ± 


0.43 


2.4 ± 1.0 


149 


13.3 


5.3 ±3.4 


1.0 ± 0.7 


J1709-4429 


r 


149.8 ± 4.1 


124 ± 2.6 


1 


.70 ± 


0.04 


4.9 ± 0.4 


16009 


373.6 


290 - 1900 


0.09 - 0.57 


J1718-3825 


r 


9.1 ± 5.8 


6.7 ± 1.9 


1 


.26 ± 


0.74 


1.3 ± 0.6 


105 


19.7 


120 ± 80 


0.09 ± 0.06 


J1732-31 


g 


25.3 ± 3.0 


24.2 ± 1.4 


1 


.27 ± 


0.14 


2.2 ± 0.3 


1002 


131.2 






J1741-2054 


g 


20.3 ± 2.0 


12.8 ± 0.8 


1 


.39 ± 


0.17 


1.2 ± 0.2 


935 


92.6 


2.2 ± 1.3 


0.24 ± 0.14 


J1744-1134 


m 


4.3 ± 1.6 


2.8 ± 0.6 


1 


.02 ± 


0.71 


0.7 ± 0.4 


78 


20.0 


0.43 ± 0.13 


0.1 ± 0.03 


J1747-2958 


r 


18.2 ± 4.2 


13.1 ± 1.7 


1 


.11 ± 


0.34 


1.0 ± 0.2 


213 


59.3 


63 - 390 


0.02 - 0.16 


J1809-2332 


g 


49.5 ± 3.0 


41.3 ± 1.6 


1 


.52 ± 


0.07 


2.9 ± 0.3 


3451 


201.9 


140 ± 140 


0.33 ± 0.33 


J1813-1246 


g 


28.1 ± 3.5 


16.9 ± 1.3 


1 


.83 ± 


0.14 


2.9 ± 0.8 


482 


39.7 






J1826-1256 


g 


41.8 ± 4.1 


33.4 ± 1.8 


1 


.49 ± 


0.11 


2.4 ± 0.3 


1152 


138 







Table 4 — Continued 



PSR 


Type a 


Fnoton I 1 lux (rioo) 
(10" 8 ph cm- 2 s- 1 ) 


Energy Mux (Giooj 
(10 -11 erg cm" 2 s" 1 ) 


r 






^-cutoff 

(GeV) 


TS 


TScutoft 


Luminosity 
(10 33 ergs" 1 ) 


hjlTicicncy 
(fn = 1) 


J1833-1034 


r 


20.5 ± 4.6 


10.1 ± 1.4 


2.24 ± 





,18 


7.7 ± 4.8 


110 


4.9 


270 ± 60 


0.01 ± 0.002 


J1836+5925 d 


g 


65.6 ± 1.8 


59.9 ± 1.3 


1.35 ± 





.04 


2.3 ± 0.1 


20982 


674.6 


<46 


<4.0 


J 1907+06 


g 


40.25 ± 3.8 


27.5 ± 1.6 


1.84 ± 





10 


4.6 ± 1.0 


1209 


59.3 






J1952+3252 


r 


17.6 ± 1.9 


13.4 ± 0.9 


1.75 ± 





12 


4.5 ± 1.2 


1008 


36.4 


64 ± 32 


0.02 ± 0.01 


J1958+2846 


g 


7.65 ± 1.6 


8.45 ± 0.83 


0.77 ± 





.31 


1.2 ± 0.2 


491 


89.2 






J2021+3651 


r 


67.35 ± 4.4 


47.0 ± 1.8 


1.65 ± 





.07 


2.6 ± 0.3 


3138 


223.5 


250 + 500 
^OU _240 


07 + ' 15 
u - u ' -0.07 


J2021+4026 d 


g 


152.6 ± 4.9 


97.6 ± 2.0 


1.79 ± 





.04 


3.0 ± 0.2 


10180 


331.4 


260 ± 150 


2.2 ± 1.3 


J2032+4127 


g 


6 ± 2.3 


11.1 ± 1.4 


0.68 ± 





.46 


2.1 ± 0.6 


487 


56.3 


34 - 170 


0.13 - 0.64 


J2043+2740 


r 


2.41 ± 0.90 


1.55 ± 0.32 


1.07 ± 





.66 


0.8 ± 0.3 


79 


15.1 


6.0 ± 3.8 


0.09 ± 0.06 


J2124-3358 
J2229+6114 


in 
r 


1.95 ± 0.49 
32.6 ± 2.2 


2.75 ± 0.42 
22.0 ± 1.0 


1.05 ± 
1.74 ± 






.34 
.08 


2.7 ± 1.0 
3.0 ± 0.5 


226 
1929 


22.9 
96.0 


21 +0 ' 42 

u - z± -0.14 

17 - 1100 


05 + 11 
0.001 - 0.05 


J2238+59 


g 


6.8 ± 1.8 


5.44 ± 0.71 


1.00 ± 





.43 


1.0 ± 0.3 


219 


37.2 







"Types are r=radio-selected, g=gamma-ray-selected, m=millisccond. 
''Here, /q is assumed to be 1, which can result in an efficiency > 1. 
c For the Crab the spectral parameters come from lAbdo et al 

d For J1836+5925 and J2021+4026 the spectral parameters come from the phase-averaged analysis (see Section 2.2). 



Note. — Results of the unbinncd maximum likelihood spectral fits for the LAT gamma-ray pulsars (see Section 2.2). Columns 3 and 4 list the on-pulse photon 
flux -Fioo and on-pulse energy flux Gioo respectively. The fits used an exponentially cutoff power-law model (see Eq. 4) with photon index T and cutoff energy 
^cutoff given in columns 5 and 6. The systematic uncertainties on -Fioo, GioO; and T due to uncertainties in the Galactic diffuse emission model have been added 
in quadrature with the statistical errors. Uncertainties in the instrument response induce additional biases of <5i<\oo = (+30%, —10%), <5Gioo = (+20%, —10%), 
5F = (+0.3, —0.1), and 8E cuto ff = (+20%, —10%). The test statistic (TS) for the source significance is provided in column 7. The significance of an exponential 
cutoff (as compared to a simple power-law) is indicated by T£ cuto ff in column 8, where a value < 10 indicates that the two models are comparable. The total 
gamma-ray luminosity L 7 and the resulting calculated gamma-ray conversion efficiency r]^ = L^/E (where fa = 1 as described in Section 3.2) are listed in 
columns 9 and 10, respectively. The uncertainties in L 7 and r\ include the flux and distance uncertainties. Nevertheless, the strong dependence of these variables 
on the measured distance (see Table [5j and beaming factor means that they should be considered with care. 
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Table 5. Pulsar distance estimates 



Pulsar Name 


Distance (kpc) 


Method 11 


(Ref 6 ) 


J0007+7303 


1.4±0.3 


K 


(30) 


J0030+0451 


0.300±0.090 


P 


(23) 


J0205+6449 


2.6-3.2 


K 


(14,32) 


J0218+4232 


2.5-4 





(1) 


J0248+6021 


2-9 





(6) 


J0437-4715 


0.1563±0.0013 


P 


(9) 


J0534+2200 


2.0±0.5 





(35) 


J0613-0200 


0.481°;!? 


P 


(18) 


J0631+1036 


0.75-3.62 





(39) 


J0633+1746 


250+ - 120 
u.zou_ 062 


p 


(11) 


J0659+1414 


n 9SO+0 033 

u.z»»_ 027 


p 


(2) 


J0742-2822 


z - u '-1.07 


DM 


(33) 


J0751+1807 


0-61S1 


P 


(28) 


J0835-4510 


0.2871°;°*? 


P 


(10) 


J1028-5819 


2.33±0.70 


DM 


(19) 


J1048-5832 


2.71±0.81 


DM 


(21) 


J1057-5226 


0.72 ±0.2 


DM 


(37) 


J1124-5916 


4.819J 


O 


(13) 


J1418-6058 


2-5 





(27,38) 


J1420-6048 


5.6±1.7 


DM 


(8) 


J1509-5850 


2.6±0.8 


DM 


(24) 


J1614-2230 


1.27±0.39 


DM 


(7) 


J1709-4429 


1.4-3.6 


O 


(26,33) 


J1718-3825 


3.82±1.15 


DM 


(25) 


J1741-2054 


0.38±0.11 


DM 


(3) 


J1744-1134 


U - J0 ' -0.035 


P 


(34) 


J1747-2958 


2-5 





(5,12) 


J1809-2332 


1.7±1.0 


K 


(29) 


J1833-1034 


4.7±0.4 


K 


(4) 


J1836+5925 


<0.8 





(16) 


J1952+3252 


2.0±0.5 


K 


(15) 


J2021+3651 


2 1+ 2 ' 1 





(36) 


J2021+4026 


1.5±0.45 


K 


(22) 


J2032+4127 


1.6-3.6 





(3) 


J2043+2740 


1.80±0.54 


DM 


(31) 


J2124-3358 


25+ 25 
u - zo -0.08 


P 


(18) 


J2229+6114 


0.8-6.5 





(17,20) 



a K distance evaluation from kinematic model; P 
from parallax; DM from dispersion measure using the 
ICordes fc Lazid l l2002h model; O from other measurements. 
For DM measurements, we assume a minimum distance un- 
certainty of 30%, as discussed in Section 3.1. 

'For DM, the reference gives the DM measurement. 

Note. — The best known distances of 37 pulsars detected 
by Fermi. Nine of the pulsars in the catalog have no dis- 
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tance estimate and arc not included in this table. 



Referenc es. — (l)lBassa et al] ||2003|); (2)lBrisken et al.l 
11200^'!; (31 ICamilo et al.l 1 2009h ; ( 4) ICamilo et al] 



(5) ICamilo et al.l l l2002b!^ (6) 



(7) ICrawford et al.l bOOrj l; (8 



(9) iDeller et al 



Fahcrl^^^al 



]200a i; (10^ 
2007h: (12) iGaensler et alJ 



CognaroN^al 



D'Ami co et al 
>odson et al. | 1 120031) 
20041 ) 



Gonzalez fc Safi-Harbl j2003l1; (14) iGreen fc Gul 
(15) lGreidanus fc Strornl ]l990l ); (16) lHalpei 



(17) lHalpern et al 



( 19) iKeith et al 
Johnstoii*^t*^n~lll99 



18 



Hot an et al 



.. ll2001al) 

(2008'); (2 0) iKothes et al 



Manchester et al. 
Ng et al] l|2q05l); 



2001 



1281 



(21) 
(23) 
(25) 

); (26)lMcGowan et al.l | |2004|); (27) 
Nice et al. | ||2005|): (29) lOka et al] 



(22) lLandecker et al] l(l98C 
(24) iManchester et al] 



19991); (30) IPineault etaL J 19931) ; ( 31)lRav et al.l | l99' 



(32) 



Roberts et al.l 




(34) [ Toscano et al 

I Van Etten et al] 1I2OO8I); (37) IWeltevrede 



(1993); 
Jl999|); 



(33) iTavlor et al 
(35) iTrimblj 



(I197S ) ; (36 ) 
Wright] j2009l); 



38) ^^^ch^aroglu^^Roman il l|l997l ): (39) IZepka' et al] 
jl996l ) 
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Table 6. Positional associations with known GeV and TeV sources for LAT-detected pulsars 



PSR 



Alt. name LAT BSL association" 



GeV associations 



Other associations 



J0007+7303 



J0030+0451 
J0205+6449 

J0218+4232 
J0248+6021 
J0357+32 
J0437-4715 

J0534+2200 Crab 

PSR B0531+21 



J0613-0200 
J0631+1036 
J0633+0632 



OFGL J0007.4+7303 
OFGL J0030.3+0450 

OFGL J0357.5+3205 
OFGL J0534.6+2201 



OFGL J0613.9-0202 
OFGL J0631. 8+1034 
OFGL J0633.5+0634 



J0633+1746 Gcminga OFGL J0634.0+1745 



J0659+1414 PSR B0656+14 
J0742-2822 PSR B0740-28 
J0751+1807 

J0835-4510 Vela OFGL J0835.4-4510 

PSR B0833-45 



J1028-5819 ■ ■ ■ OFGL J1028.6-5817 

J1048-5832 PSR B1046-58 OFGL J1047.6-5834 

J1057-5226 PSR B1055-52 OFGL J1058.1-5225 
J1124-5916 

J1418-6058 ■ ■ ■ OFGL J1418.8-6058 

J1420-6048 



3EG J0010+7309 
EGR J0008+7308 
GEV J0008+7304 
1AGL J0006+7311 
EGR J0028+0457 



3EG J0534+2200 
EGR J0534+2159 
GEV J0534+2159 
1AGL J0535+2205 



3EG J0631+0642 
EGR J0633+0646 
GEV J0633+0645 
3EG J0633+1751 
EGR J0633+1750 
GEV J0634+1746 
1AGL J0634+1748 



3EG J0834-4511 
EGR J0834-4512 
GEV J0835-4512 
1AGL J0835-4509 
3EG J1027-5817 
GEV J1025-5809 
3EG J1048-5840 
EGR J1048-5839 
GEV J1047-5840 
3EG J1058-5234 
EGR J1058-5221 
GEV J1059-5218 
1AGL J1058-5239 



3EG J1420-6038 
GEV J1417-6100 
1AGL J1419-6055 
3EG J1420-6038 



SNR CTA 1 2 
PWN G119.5+10.2 1 



SNR/PWN 3C 58 2 
PWN G119.5+10.2 1 



PWN G253.4-42.0 1 
SNR/PWN G184.6-5.8 
HESS J0534+220 4 



PWN G195.1+4.3 1 
MGRO J0632+17 3 > n 



SNR 203.0+12.0 



SNR/PWN G263. 9-3.3 1 ' 2 
HESS J0835-455 5 



PWN G287.4+0.58 1 



MSH 11-54 

SNR/PWN G292.0+1.8 1 ' 2 
PWN G313.3+0.1 1 
HESS J1418-609 6 

PWN G313.6+0.3 1 
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Table 6 — Continued 



PSR 



Alt. name 



LAT BSL association" 



GeV associations 



Other associations 



J1459-60 
J1509-5850 c 
J1614-2230 
J1709-4429 



J1718-3825 
J1732-31 



J1741-2054 
J1744-1134 
J1747-2958 
J1809-2332 



J1813-1246 
J1826-1256 



J1833-1034 
J1836+5925 

J1907+06 



OFGL J1459.4-6056 
OFGL J1509.5-5848 

PSR B1706-44 OFGL J1709. 7-4428 



OFGL J1732.8-3135 



OFGL J1742.1-2054 



OFGL J1809.5-2331 



OFGL J1813.5-1248 
OFGL J1825.9-1256 



OFGL J1836.2+5924 



OFGL J1907.5+0602 



J1952+3252 d PSR B1951+32 OFGL J1953.2+3249 



J1958+2846 
J2021+3651 
J2021+4026 
J2032+4127 



J2043+2740 
J2124-3358 
J2229+6114 



OFGL J1958.1+2848 
OFGL J2020.8+3649 
OFGL J2021. 5+4026 
OFGL J2032.2+4122 



OFGL J2124.7-3358 
OFGL J2229.0+6114 



EGR J1418-6040 
GEV J1417-6100 
1AGL J1419-6055 

1AGL J1506-5859 
3EG J1616-2221 
3EG J1710-4439 
EGR J1710-4435 
GEV J1709-4430 
1AGL J1709-4428 

3EG J1734-3232 
EGR J1732-3126 
GEV J1732-3130 
3EG J1741-2050 

1AGL J1746-3017 
3EG J1809-2328 
GEV J1809-2327 
1AGL J1809-2333 
GEV J1814-1228 
3EG J1826-1302 
GEV J1825-1310 
1AGL J1827-1277 



3EG J1835+5918 
GEV J1835+5921 
1AGL J1836+5923 
GEV J1907+0557 
1AGL J1908+0613 



3EG J1958+2909 
GEV J1957+2859 
GEV J2020+3658 
1AGL J2021+3652 
3EG J2020+4017 
1AGL J2021+3652 
3EG J2033+4118 
EGR J2033+4117 
1AGL J2032+4102 



3EG J2227+6122 
EGR J2227+6114 



HESS J1420-607 



PWN G319.97-0.62 1 

SNR/PWN G343. 1-2.3 x ' 2 
HESS J1708-443 7 



HESS J1718-385 



PWN G359.23-0.82 1 
PWN G7.4-2.0 1 



PWN G18. 5-0.4 1 



SNR/PWN G21.5-0.9 l > 2 
HESS J1833-105 9 



MGRO J1908+063 3 
HESS J1908+06 10 
SNR CTB 80 2 
PWN G69.0+2.7 1 



PWN G75. 2+0.1 1 
MGRO J2019+37 3 
SNR 7 Cygni 2 
SNR G78.2+2.1 2 
MGRO J2031+41 3 



PWN G10.9-45.4 1 
PWN G106.6+2.9 1 
MGRO J2228+61 3 > n 



Table 6 — Continued 



PSR Alt. name LAT BSL association" GeV associations 6 Other associations 

GEV J2227+6101 
1AGL J2231+6109 

J2238+59 



"Source designator from the LAT Bright Source List l|Abdo et al.ll2009il) . 

^Source designator(s) from the 3rd EGRET (3EG: lHartman et al.ll999h. Revise d EGRET (EGR: 
ICasandiian fc Grenierll2008h High -energy EGRET (GEV lLamb fc Macomblll997h and/or the first 
AGILE C1AGL: IPittori et al.ll2009h catalogs. 

C PSR J1509-5850 should not be confused with PSR B1509-58 observed by CGRO llKuiper et al.l 

Il999h . 

d Whilc pulsations from PSR J1952+3252 were detected in EGRET data, it was never cataloged 
as a point source. 



Note. — Alternate names for the pulsars in this catalog are given in column 2. Positional 
associations with SNRs, PWNe and selected TeV sources are provided in column 5. 



References. — I- iRoberts et all d2005l), 2. iGreenl j2009h. 3. lAbdo et al.l <2009lj|. 4 
lAharonian et al.l i2006d). 5 . Aharonian et al.l i2006bl). 6. lAharonian et al. N2006af). 7. iHoppe et al.l 
1 20091). 8. Aharonian et al.l J2007tl . 9. iDiannati-Atai et al.l (120071 1 . 10. lAharonTan et ah I J2009l l. 11. 



iGoodman &: Sinnij J2009F 
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Fig. 1. — Pulsar sky map in Galactic coordinates. Blue squares: gamma-ray-selected pulsars. 
Red triangles: millisecond gamma-ray pulsars. Green circles: all other radio loud gamma-ray pul- 
sars. Black dots: Pulsars for which gamma-ray pulsation searches were conducted using rotational 
ephemerides. Gray dots: Known pulsars which were not searched for pulsations. 
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Fig. 2. — P — P diagram. Dashed lines: characteristic age r c . Dot-dashed lines: rotational energy 
loss rate E. Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray 
pulsars. Green circles: all other radio loud gamma-ray pulsars. Black dots: Pulsars for which 
gamma-ray pulsation searches were conducted using rotational ephemerides. Gray dots: Known 
pulsars which were not searched for pulsations. 
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Fig. 3. — Galactic plane pulsar distribution (polar view). The star represents the Galactic center. 
The two circles centered at the Earth's position have radii of 3 kpc and 5 kpc. For pulsars with dif- 
ferent possible distances, the nearer values from Table [5] are used. Note that the millisecond pulsars 
(MSPs) , while having a significantly lower E than the other pulsars (see Figure [8]) , are detectable 
due to their close proximity. The one exception (PSR J02 18+4232) also exhibits a significantly 
higher E than the other MSPs. Blue squares: gamma-ray-selected pulsars. Red triangles: millisec- 
ond gamma-ray pulsars. Green circles: all other radio loud gamma-ray pulsars. Black dots: Pulsars 
for which gamma-ray pulsation searches were conducted using rotational ephemerides. Gray dots: 
Known pulsars which were not searched for pulsations. 
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Fig. 4. — Phase difference A between the gamma-ray peaks, versus the phase lag 5 between the 
main radio peak and the nearest gamma-ray peak. Pulsars without a radio detection are plotted 
with 5 = 0. With present light curves we cannot generally measure A < 0.15; objects classified as 
single-peaked are plotted with A=0. Two such objects, both MSPs, are off the plot at 5 > 0.8. 
Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray pulsars. Green 
circles: all other radio loud gamma-ray pulsars. The plot has been corrected for the updated value 
of 5 for PSR J1124-5916, as per an Erratum sent to the ApJ (December 2010). 
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Fig. 5. — Magnetic field strength at the light cylinder £?lc versus pulsar characteristic age r c . 
Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray pulsars. Green 
circles: all other radio loud gamma-ray pulsars. Black dots: Pulsars for which gamma-ray pulsation 
searches were conducted using rotational ephemerides. Gray dots: Known pulsars which were not 
searched for pulsations. 
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Fig. 6. — Gamma-ray luminosity L 7 versus the rotational energy loss rate E. Dashed line: L 7 
equal to E. Dot-dashed line: L 7 proportional to the square root of E. L 7 is calculated using a 
beam correction factor /q = 1 for all pulsars and the integral energy flux Gioo from the on-pulse 
spectral analysis (see Section 2.2), except for PSRs J1836+5925 and J2021+4026 which use the total 
background-corrected phase- averaged flux, including a relatively bright unpulsed component (see 
Section 2.2). For the Crab we also plot the total high energy luminosity, L tot = Lx +L~, indicated 
by *. Several notable pulsars have been labeled. Blue squares: gamma-ray-selected pulsars. Red 
triangles: millisecond gamma-ray pulsars. Green circles: all other radio loud gamma-ray pulsars. 
Unfilled markers indicate pulsars for which only a DM-based distance estimate is available (see 
Table [5]). Pulsars with two distance estimates have two markers connected with dashed error bars. 
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Fig. 7. — Value of the exponential cutoff E cuto g versus the magnetic field at the light cylinder, -Blc- 
The statistical uncertainties on E cuto Q are shown. An additional systematic bias of (+20%, —10%) 
may affect E cuto fi (see text). The histogram of -^cutoff values is projected along the right-hand axis. 
Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray pulsars. Green 
circles: all other radio loud gamma-ray pulsars. 
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Fig. 8. — Photon index T versus the rotational energy loss rate, E. For T, the statistical uncertain- 
ties combined with the systematic uncertainties due to the diffuse emission model are shown. An 
additional systematic bias of (+0.3, —0.1) affects T (see text). The histogram of the photon indices 
is projected along the right-hand axis. Blue squares: gamma-ray-selected pulsars. Red triangles: 
millisecond gamma-ray pulsars. Green circles: all other radio loud gamma-ray pulsars. 
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Fig. 9. — Aitoff projection sky map of the 5a sensitivity in units of logarithmic photon flux 
(Log(L 7 ) ph cm" 2 s^ 1 ) for six months of Fermi LAT sky-survey data. The sensitivity analysis 
uses the model of the diffuse gamma-ray background described in the text (Section 4), and pulsar 
spectra with differential photon indices of T = 1.4 with an exponential cutoff energy of E cuto f[ = 2.2 
GeV. 




Fig. 10. — Measured integral photon flux above 100 MeV, -Fioo, versus the 5<r flux sensitivity de- 
scribed in Figure [9l For Fioq, the statistical uncertainties combined with the systematic uncertain- 
ties due to the diffuse emission model are shown. An additional systematic bias of (+30%, —10%) 
affects -Fioo (see text). The effective blind search sensitivity is comparable to the 2 x 5cr line, 
although a few pulsars are discovered at lower flux, presumably due to favorable pulse profiles, 
spectra or local backgrounds. Blue squares: gamma-ray-selected pulsars. Red triangles: millisec- 
ond gamma-ray pulsars. Green circles: all other radio loud gamma-ray pulsars. 
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Fig. 11. — LogN-LogS distribution as described in Section 4 for all the detected pulsars (black 
dashed line), the radio-selected gamma-ray pulsars including MSPs (grey histogram), and the 
gamma-ray-selected pulsars (blue hatched histogram). 
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Fig. 12. — Pulsar "detectability" metric (E 1 / 2 /d 2 , normalized to Vela) vs. spin period. Detected 
MSPs (red triangles) and young pulsars (radio-selected, green circles; gamma-ray-selected, blue 
squares) all have high values of this metric. For objects with a distance range in Table [5j we use 
here the geometric mean of the maximum and minimum values. Searched, but presently undetected 
objects (gray dots) are plotted using DM-derived distances. For the possible causes of non-detection 
see Section 5.1. 
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Fig. 13. — Separations A between the gamma-ray peaks, for those pulsars with two identified 
peaks, versus the spin-down power E. The histogram of peak separations is projected along the 
right-hand axis. Blue squares: gamma-ray-selected pulsars. Red triangles: millisecond gamma-ray 
pulsars. Green circles: all other radio loud gamma-ray pulsars. 
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Appendix: Gamma-ray Pulsar Light Curves 




Fig. A-l.— Light curves for PSR J0007+7303 (P 
316 ms). 
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Fig. A-2. 
4.87 ms). 



Light curves for PSR J0030+0451 (P 




Fig. A-3.— Light curves for PSR J0205+6449 (P 
65.7ms). 




Fig. A-4— Light curves for PSR J0218+4232 (P 
2.32 ms). 




Fig. A-5.— Light curves for PSR J0248+6021 (P 
217ms). 




A-6.— Light curves for PSR J0357+32 (P = 444 




Fig. A-7— Light curves for PSR J0437-4715 (P 
5.76 ms). 




Fig. A-8.— Light curves for PSR J0534+2200 (P = 
33.1ms, Crab pulsar). The zero of phase is set to the 
radio precursor. 




Fig. A-9.— Light curves for PSR J0613-0200 
3.06 ms). 
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Fig. A-10.— Light curves for PSR J0631+1036 (P 
288 ms). 




Fig. A- 11.— Light curves for PSR J0633+0632 (P 
297 ms). 




Fig. A-12.— Light curves for PSR J0633+1746 (P 
237ms, Geminga pulsar). 




Fig. A-13.— Light curves for PSR J0659+1414 (P 
385 ms). 




Fig. A-14.— Light curves for PSR J0742-2822 (P 
167 ms). 





Fig. A-15.— Light curves for PSR J0751+1807 (P 
3.48 ms). 




Fig. A-16.— Light curves for PSR J0835-4510 (P 
89.3 ms, Vela pulsar). 




Fig. A-17— Light curves for PSR J1028-5819 (P 
91.4 ms). 




Fig. A-18.— Light curves for PSR J1048-5832 (P 
124 ms). 




Fig. A-19 — Light curves for PSR J1057-5226 (P 
197ms, PSR B1055-52). 




Fig. A-20.— Light curves for PSR Jl 124-5916 (P = 
135 ms). The radio phase alignment has been corrected, 
as per an Erratum sent to the ApJ (December 2010). 




Fig. A-21.— Light curves for PSR J1418-6058 (P 
111 ms). 




Fig. A-22.— Light curves for PSR J1420-6048 (P 
68.2 ms). 




Fig. A-23.— Light curves for PSR J1459-60 (P 
103 ms). 




Fig. A-24.— Light curves for PSR J1509-5850 (P 
88.9 ms). 
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Fig. A-25.— Light curves for PSR J1614-2230 (P 
3.15 ms). 




Fig. A-26.— Light curves for PSR J1709-4429 (P 
102 ms, PSR B1706-44). 
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Fig. A-27. 
74.7 ms). 



Light curves for PSR J1718-3825 (P 




Fig. A-28.— Light curves for PSR J1732-31 
197 ms). 
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Fig. A-29.— Light curves for PSR J1741-2054 (P 
414 ms). While t his pulsar is detected in the radio 



(jCamilo et al.l 12009 ). it was discovered by the LAT and 
is considered a gamma-ray-selected pulsar. 




Fig. A-30.— Light curves for PSR J1744-1134 (P 
4.08 ms). 




Fig. A-31.— Light curves for PSR J1747-2958 (P 
98.8 ms). 




Fig. A-32.— Light curves for PSR J1809-2332 (P 
147 ms). 
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Fig. A-33.— Light curves for PSR J1813-1246 (P 
48.1 ms). 
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Fig. A-34.— Light curves for PSR J1826-1256 (P 
110 ms). 




Fig. A-35.— Light curves for PSR J1833-1034 (P 
61.9 ms). 




Fig. A-36.— Light curves for PSR J1836+5925 (P 
173 ms). 




Fig. A-37.— Light curves for PSR J1907+06 (P 
107 ms). 
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Fig. A-38.— Light curves for PSR J1952+3252 (P 
39.5 ms, PSR B1951+32). 




Fig. A-39.— Light curves for PSR J1958+2846 (P 
290 ms). 




Fig. A-40.— Light curves for PSR J2021+3651 (P 
104 ms). 




Fig. A-41 .— Light curves for PSR J2021+4026 (P 
265 ms) . 
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Fig. A-42.~ Light curves for PSR J2032+4127 (P = 
14 3 ms). While t his pulsar is detected in the radio 



(jCamilo et al.l 12009 ). it was discovered by the LAT and 
is considered a gamma-ray-selected pulsar. 




Fig. A-43.— Light curves for PSR J2043+2740 (P 
96.1 ms). 





Fig. A-45.— Light curves for PSR J2229+6114 (P 
51.6 ms). 




Fig. A-46.— Light curves for PSR J2238+59 (P 
163 ms). 



